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Abstract 

The  advent  of  SD-periodic  dielectric  structures  possessing  photonic  band 
gaps  (PBGs)  and  the  ability  to  incorporate  defects  in  the  structure,  wherein 
localized  high-Q  modes  are  formed  at  PBG  frequencies,  should  offer  far  greater 
control  over  (i)  the  spatial  modulation  of  the  field  amplitude,  and  (ii)  the  spec¬ 
tral  distribution  of  its  mode  density.  We  have  been  exploring  the  following 
novel  quantum  optical  processes  which  rely  on  the  ’’design”  of  the  afore¬ 
mentioned  field  characteristics  in  PBG  structures:  (a)  Quantum  states  via 
nonadiabatic  periodic  transitions:  Unusual  nonadiabatic  dynamics  of  field- 
dressed  atomic  states  is  obtained  for  atoms  moving  along  a  near-resonant 
periodically-modulated  field  mode.  Such  dynamics,  followed  by  measure¬ 
ment  of  the  atomic  excitation,  permits  an  extremely  powerful  control  over 
the  photon-number  and  phase  distributions  of  the  field,  (b)  Pump-mode 
Fock-state  generation:  Atomic  states  dressed  by  the  defect  field  can  decay 
via  coupling  to  the  mode  continuum  into  a  state,  which  upon  measurement 
of  the  atomic  state,  corresponds  to  a  single  Fock  state  of  the  field,  (cl)  Las¬ 
ing  without  inversion  -  atomic  coherence  by  spontaneous  decay:  Spontaneous 
decay  of  an  atom  whose  resonance  is  located  between  two  PBGs  is  found  to 
produce  coherence  between  adjacent  atomic  sublevels.  This  coherence  can 
be  used  to  generate  electromagnetically-induced  transparency,  lasing  without 
inversion  and  nonclassical  light.  (c2)  Two-atom  interactions  and  coherence: 
We  predict  drastic  suppression  or  enhancement  of  the  resonant  excitation 
transfer  and  the  cooperative  frequency  shift  in  systems  of  two  identical  atoms 
placed  within  PBG  structures,  both  at  near-zone  and  far- zone  inter-atom 
separations.  These  effects  stem  from  the  competition  of  PBG-induced  single¬ 
atom  spectral  shifts  and  the  inter-atom  resonant  dipole-dipole  interaction. 
The  resulting  two-atom  stable  correlated  states  can  be  used  for  lasing  with¬ 
out  inversion,  (d)  Near-resonant  gap  solitons  and  optical  ^excitons”:  We 
demonstrate  that  pulse  transmission  through  near-resonant  media  embedded 
in  ID-periodic  dielectric  structures  can  yield  a  novel  type  of  gap  solitons, 
which  allows  for  near-resonant  seif-induced  transparency  at  PBG  frequen¬ 
cies.  At  the  QED  level,  we  predict  the  existence  of  multiphoton  bound  states, 
which  are  the  optical  analogs  of  multi-excitons  at  PBGs  of  nonlinear  periodic 
structures. 


Keywords:  Photonic  band  structures,  quantum  optics,  quantum  state  preparation,  lasing 
without  inversion,  cooperative  (collective)  effects,  coherence  effects,  gap  solitons. 


I.  QUANTUM  FIELD  STATES  VIA  NON- ADIABATIC  TRANSITIONS 


A.  Background 

PBG  structures  are  expected  to  allow  better  control  over  one  of  the  basic  properties 
that  determine  the  field-atom  interaction,  i.e.,  the  spatial  volume  and  shape  of  narrow- 
linewidth  (high-Q)  modes.  Many  of  the  investigated  structures  are  composed  of  void  spaces 
periodically  alternating  with  dielectric  material.  Atomic  beams  collimated  to  a  micron- 
size  waist  can  therefore  pass  through  void  channels  in  structures  possessing  PBGs  at  near- 
infrared  wavelengths.  The  resulting  field-atom  interactions  are  expected  to  extend  the  scope 
of  QED  processes  realizable  in  the  regime  of  a  single  high-Q  mode  near  resonance  with 
atomic  transition,  reversible  energy  exchange  is  obtainable  between  the  atom  and  field.  In 
the  simple  case  of  an  atom  traveling  in  a  spatially  uniform  high-Q  mode,  this  exchange 
is  described  by  the  fundamental  Jaynes-Cummings  model  (JCM)  [1,2].  Atomic  motion 
through  a  periodically  modulated  high-Q  mode,  in  a  Fabry-Perot  resonator  or  a  defect  at  a 
PBG  frequency,  has  been  predicted  to  drastically  modify  the  resonant  JCM  dynamics,  and, 
as  shown  by  us,  allow  better  control  over  the  generation  of  sub-Poissonian  photon  statistics 
in  such  a  mode  [3].  In  the  present  study,  we  have  posed  the  question:  can  photon  statistics 
and  other  nonclassical  properties  of  the  field  be  better  controlled  by  near-resonant  periodic 
modulation  of  the  field-atom  coupling? 


B.  Results 

A  novel  mechanism  for  generating  nonclassical  field  states  has  been  discovered  (item  6  - 
Sec.V).  based  on  non- adiabatic  periodic  modulation  of  the  parameters  governing  the  field- 
atom  coupling.  This  mechanism  has  been  investigated  in  the  context  of  a  two-level  atom 
moving  along  a  spatially-periodic  near- resonant  field  mode  in  a  PBG  structure,  resulting  in 
temporal  periodicity  of  the  Rabi  frequency.  We  have  shown  that  steady  growth  with  number 
of  periods  is  possible  for  the  field-dressed  atomic  transition  probability.  Such  steady  growth 
(parametric  amplification)  is  obtained  whenever  the  probability  amplitudes  in  a  single  period 
interfere  constructively,  which  corresponds  to  the  adiabatic  modulation  phase  in  each  period 
being  an  integer  multiple  of  tt. 

An  advantageous  property  of  such  system  is  the  possibility  of  selective  generation  of  spe¬ 
cific  nonclassical  states  of  the  field:  (a)  For  an  atom  interacting  with  a  spatially  sinusoidal 
quantized  field,  the  possibility  of  strong  enhancement  of  the  transition  probability  (by  con¬ 
structive  interference  at  many  nodes  of  the  field)  is  very  sensitive  to  the  number  of  photons 
in  the  quantized  field.  Hence,  the  transition  amplitudes  may  interfere  constructively  only 
for  certain  photon-number  states.  Thus  a  successful  measurement  of  the  atomic  excitation 
results  in  nearly  a  Fock  state  of  the  field,  if  only  one  photon  number  satisfies  the  construc¬ 
tive  interference  condition  for  the  transition.  The  regime  in  which  several  number  states 
are  singled  out  from  the  initial  distribution  by  the  constructive  interference  and  the  atomic- 
excitation  measurement  results  in  the  generation  of  amplitude  or  phase  squeezed  states,  (b) 
The  periodic  modulation  of  the  field-atom  detuning  controlled  by  sufficiently  many  param¬ 
eters  gives  us  the  possibility  to  generate  an  arbitrary  preselected  superposition  of  photonic 
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number  states  (quantum  st,ate  engineering).  The  number  of  components  in  the  preselected 
superposition  is  determined  by  the  number  of  control  parameters  (Appendix  A). 


II.  PUMP-MODE  FOCK-STATE  GENERATION  BY  RESONANCE 
FLUORESCENCE  IN  PHOTONIC  BAND  STRUCTURES 

A.  Background 

Resonance  fluorescence  (RF)  is  the  emission  or  scattering  of  photons  into  a  mode  contin¬ 
uum  by  a  two-level  atom  which  is  driven  by  a  near-resonant  single-mode  field  [4].  This  has 
been  one  of  the  most  extensively  studied  processes  in  quantum  optics,  through  which  several 
fundamental  features  of  field-atom  interaction  have  been  revealed,  such  as  the  Mollow  triplet 
spectrum  [5]  and  nonclassical  field  correlations,  manifested  by  antibunching  [6],  quadrature 
squeezing  [7]  and  sub-Poissonian  photon  statistics  [8].  More  recently,  the  interest  in  RF  has 
shifted  to  its  modifications,  taking  place  when  the  mode  continuum,  into  which  the  driven 
atom  radiates, differs  from  ordinary  free-space  vacuum.  Thus,  a  squeezed  ’’vacuum”  contin¬ 
uum,  in  which  one  field  quadrature  of  pairwise-correlated  modes  exhibits  reduced  quantum 
fluctuations,  has  first  been  shown  to  affect  spontaneous  decay  [9],  and  subsequently  has  been 
demonstrated  to  modify  the  Mollow-triplet  peaks,  depending  on  the  squeezing  bandwidth 
("coloring”  of  fluctuations)  and  the  phase  between  the  squeezed  vacuum  and  coherent  driv¬ 
ing  fields  [10].  Analogous  modifications  of  RF  have  been  revealed  in  resonator  structures 
with  ’’coloring”  of  the  continuum  mode  density  [11].  These  modifications  are  expected  to  be 
especially  dramatic  in  periodic  dielectric  structures  that  exhibit  photonic  band  gaps  (PBGs), 
i.e.,  forbidden  spectral  bands  at  all  emission  angles  [12].  For  a  nearly-classical  non-depleted 
driving  field,  the  inhibition  of  one  of  the  Mollow  triplet  spectral  sidebands  by  a  PBG  has 
been  shown  to  inhibit  the  other  sideband,  thereby  causng  the  disappearance  of  antibunching 
[11].  The  latter  phenomenon  has  been  analyzed  using  the  common  assumption  of  locally- 
Markovian  response  of  the  mode  continuum.  On  the  other  hand,  studies  of  spontaneous 
decay  into  strongly  colored  continua  [13]  indicate  that  non-Markovian  features  may  become 
prominent  near  an  abrupt,  nonanalytic  cutoff  of  the  continuum,  in  a  resonator  or  a  periodic 
dielectric  structure.  One  such  feature  is  the  formation  of  a  non-decaying  photon-dressed 
atomic  state  in  the  forbidden  band  [14-16]. 


B.  Results 

In  recent  months  [17]  (Item  1  -  Sec.V)  we  have  studied  the  dynamical  effect  of  RF  on  the 
photon  statistics  and  correlations  of  the  driving  field,  which  undergoes  depletion  via  dressed- 
state  decay  near  abrupt  PBG  edges.  The  most  striking  outcome  of  this  process  is  that  an 
initially  broad  (Poissonian  or  thermal)  photon-number  distribution  of  the  driving  field  can 
become  highly  sub-Poissonian,  with  a  narrow  peak  about  a  single  Fock  state,  following  a 
measurement  of  the  atomic  excitation.  The  realization  of  such  distributions  may  therefore 
be  a  major  stepping  stone  towards  the  construction  of  novel  quantum  light  sources,  including 
photon-number  states.  We  also  show  that  the  dressed-state  non-diagonal  matrix  elements 
yield,  upon  measuring  the  atomic  state,  coherences  between  Fock  states  of  the  field.  Hence, 


3 


the  two  predominant  ’’attractor”  states,  say,  |no)  and  |no  —  1),  whenever  the  defect  mode 
is  shifted  from  the  PBG  center,  will  be  partly  correlated,  if  an  initial  coherence  exists  either 
in  the  atom  or  in  the  field.  Therefore  deviations  from  phase  randomness  of  the  field  are 
anticipated. 

Another  striking  effect  of  the  stable  (non-decaying)  dressed-state  coherences  is  the  ap¬ 
pearance  of  beats  ( oscillations)  of  Fock-state  populations  of  the  field  prepared  by  measuring 
the  final  atomic  state.  These  oscillations  will  be  strongly  pronounced  in  the  case  when  the 
initial  photon  distribution  spans  many  PBGs  in  the  spectral  density  of  modes  (DOM).  In 
this  case  the  field  distribution  will  consist  of  many  partly-correlated  Fock  states  [17]. 


III.  LASING  WITHOUT  INVERSION: 

III.l  Single-Atom  Systems 
A.  Background 

If  the  lower  state  of  the  resonant  atomic  transition  consists  of  two  sublevels  |^)  and  \g') 
a  A-configuration),  then  the  atomic  state  coherence  can  be  used  in  a  scheme  of  lasing 
without  inversion(LWI)  between  the  |y)  —  j^f')  doublet  and  a  higher  atomic  level  |ri).  It  has 
been  long  known  that,  if  the  single-photon  Rabi  frequency  of  a  cavity  mode,  Oi,  is  much 
larger  than  its  linewidth  Pi,  then  the  dressed  atomic  state  will  be  split  after  time  ~ 
into  a  nearly-equal  superposition  of  two  dressed  states  whose  energies  wj  and  pushed 

in  opposite  directions  from  the  cavity-mode  peak  by  roughly  ±Ui/2.  We  have  posed  the 
question:  can  we  base  LWI  schemes  on  similar,  spontaneously  induced  coherence  in  PBG 
structures?  (Item  7  -  Sec.V). 


B.  Results:  Spontaneous  Coherence  Formation  in  PBG  Structures 

Assume  that  the  transition  between  a  ground  \g)  and  an  excited  |e)  levels  of  the  atom 
is  near  resonant  with  a  narrow  photonic  band.  The  band  is  assumed  to  be  well  separated 
from  other  bands,  so  that  the  influence  of  the  latter  can  be  neglected.  Moreover,  the  band 
is  assumed  to  be  sufficiently  dense  and  narrow,  i.e., 

>  Aft,  |A™|,  (1) 

where  Gm  and  Aft  are  the  maximum  value  and  the  width  of  the  function  G{u>)  characterizing 
the  stregth  of  the  atom-band  coupling,  A^  =  Wa  —  is  the  center  of  gravity  of  the 

band. 


1 

~~7  d(jJuG{u)^  (2) 

/i  J Uljj 

(jJu  and  are  the  lower  and  upper  band  edges  respectively,  and 

r^L 

A~  dujG{ij),  (3) 

Juu 
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Then  the  level  e  can  be  shown  [16]  to  split  to  two  levels  1  and  2,  whose  separation  to  first 
approximation  is 


E2  —  El  —  hyj 4A  +  —  hflR,  (4) 

where  Hr  is  the  vacuum  Rabi  frequency.  An  excited  atom  injected  at  t  =  0  into  the  PBS 
in  the  case  (1)  undergoes  an  insignificantly  small  spontaneous  decay  over  a  time  of  order 
l/flRi  after  which  the  excited-state  wavefunction  becomes  a  coherent  superposition  of  the 
states  I'i/’i)  and  1^2)  [16], 

\^{t))  =  v/driV^i)e-‘“>'  -f-  (5) 

where  Uj  =  Ejjh.  If  there  is  a  nonzero  matrix  element  of  the  dipole  moment,  dze-,  for  the 
transition  between  the  level  e  and  another  level  3,  then  the  dipole  moment  matrix  elements 
for  the  transitions  1-3  and  2-3  are  given  by 


dzi  —  dz2  —  V^^Se- 


(6) 


For  the  near-resonant  case,  <C  4A,  which  is  of  special  interest  here,  one  gets  in  the  first 
approximation 


1  Ula-\-  OJ-m  J-r 

ci=C2  =  -,  u;i,2  = - 2 - Tv  A.  (() 

The  density  matrix  elements  of  the  system  consisiting  of  the  levels  1  and  2  are  [cf.  Eqs.  (5) 
and  (7)] 

Y  g2iv74? 

P\\{t)  =  p22{t)  =  Pl2{t)  =  p2i{t)  =  — (8) 


C.  Results:  Lasing  Without  Inversion 

The  coherence  between  the  levels  1  and  2  can  be  utilized  to  produce  lasing  without 
inversion  [18].  Following  Refs.  [19-21],  consider  a  A-system.  The  field 

E{t)  =  Eoe-'"*  +  c.c.  (9) 

couples  the  upper  level  3  with  the  lower  levels  1  and  2.  The  levels  1,  2,  and  3  have  have 
the  reciprocal  lifetimes  respectively  71,  72,  and  73  =  73^  -1-  732  +  73,  where  73^,  732,  and  73 
are  the  spontaneous  rates  of  transitions  from  the  level  3  to  the  levels  1,  2,  and  all  others 
respectively,  whereas  the  reciprocal  lifetimes  of  the  coherences  pi2  and  P13  are  721  and  731 
respectively.  The  levels  1  and  2  can  have  finite  lifetimes  due  to  one  or  several  of  the  following 
factors:  spontaneous  transitions  to  levels  other  than  the  ground  state,  finite  time  of  presence 
of  the  atom  in  the  structure,  finite  density  of  modes  outside  the  narrow  band  (which  causes 
the  splitting  of  the  level  e),  and  collision-induced  transitions.  The  open  3- level  sysytem  is 
described  here  by  an  unnormalized  density  matrix  R{t),  whose  diagonal  elements  yield  the 
number  of  atoms  in  the  laser  cavity  in  the  respective  level. 

Our  results  show  that  a  spontaneously-formed  coherent  superposition  of  the  states  1 
and  2  allows  only  a  very  small  fraction  of  atoms  in  these  states  to  absorb  the  field  E{t). 
Therefore  there  can  be  gain  whenever  the  number  of  the  atoms  at  the  level  3  is  much  less 
than  the  number  of  the  atoms  at  the  lower  levels  1  and  2  and  iui2  ^  72i)73i  (Appendix  B). 
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III.2  Two- Atom  Effects 


A.  Background 

Studies  of  cooperative  atomic  effects  in  resonantors  have  been  extensions  of  the  Tavis- 
Cummings  model  [22].  This  model  assumes  many  atoms  identically  coupled  to  a  single  mode 
and  ignores  the  symmetry- breaking  dipole-dipole  effects,  which  are  important  at  near-zone 
separations  [23,24]. 

Our  study  is  aimed  at  gaining  new  insight  into  cooperative  effects:  a  pair  of  identi¬ 
cal  two-level  atoms  (or  excitons),  sharing  a  photon  with  one  or  many  high-Q  modes  in  a 
photonic-band  structure  or  resonator.  Such  effects  may  be  important  for  the  spontaneous 
formation  of  coherence  in  two-atom  systems  and  its  utilization  for  lasing  without  inversion 
(see  above).  We  have  developed  a  formalism  capable  of  treating  two  atoms  coupled  to  a  field 
with  an  arbitrary  xnode-density  spectrum  in  a  non-perturbative  fashion.  Such  a  formalism  is 
necessary  because  the  standard  perturbative  treatment  of  two-atom  coupling,  to  second  or¬ 
der  in  the  field  [25],  is  inadequate  for  a  reservoir  whose  mode-density  spectrum  does  not  vary 
smoothly,  as  demonstrated  already  by  our  comprehensive  theory  of  a  single  atom  coupled 
to  such  a  reservoir  [16]. 


B.  Results 

Our  treatment  (Item  4  Sec.V)  demonstrates  that  two  identical  atoms  interacting  via 
a  near-resonant  narrow-linewidth  mode  (or  degenerate  band)  and  an  off-resonant  reservoir 
can  exhibit  a  much  richer  variety  of  spectral  and  dynamical  features  than  what  is  currently 
known.  It  must  be  viewed  as  a  system  of  three  mutually  coupled  excited  states  (as  opposed 
to  two  such  states  in  previously  studied  models  [25,22]).  The  dipole-dipole  (RDDI)  coupling 
and  the  mode-induced  Rabi  splitting  are  inseparable  in  this  system.  The  interatom  coupling 
results  from  mixing  of  all  three  states,  which  is  most  striking  for  separations  where  the  two 
predominantly-populated  levels  nearly  cross.  Then  the  competing  RDDI  and  Rabi  splittings 
cause  strong  interference  of  the  symmetric  and  antisymmetric  two-atom  excited  states,  lead¬ 
ing  to  decoupling  of  single-atom  excited  states.  This  occurs  at  near-zone  (quasimolecular) 
separations  (as  small  as  10“^  A^)  and  corresponds  to  suppression  of  interatom  excitation 
transfer,  compared  to  the  electrostatic  near-zone  limit  of  RDDI. 

The  present  predictions,  particularly  the  suppression  of  excitation  transfer  at  quasimolec¬ 
ular  separations,  may  be  important  in  various  systems  of  two  atoms,  molecules,  or  excitons 
within  high-(5  Bragg  resonators  (Appendix  C). 
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,IV.  PROPAGATION  EFFECTS 

IV. 1  Self-Induced  Transparency  in  Photonic  Band  Structures;  Gap  Solitons  Near 

Absorption  Resonances 

A.  Background 

Pulse  propagation  in  a  non-uniform  resonant  medium,  e.g.,  a  periodic  array  of  resonant 
films,  can  destroy  self-induced  transparency  (SIT)  [26,27],  because  the  pulse  area  is  then 
split  between  the  forward  and  backward  (reflected)  coupled  waves,  and  is  no  longer  conserved 
[28,29].  Should  we  then  anticipate  severely  hampered  transmission  through  a  medium  whose 
resonance  lies  in  a  reflective  spectral  domain  (photonic  band  gap)  of  a  periodically-layered 
structure  (a  Bragg  reflector)?  We  have  shown  analytically  [30]  (Item  2  -  Sec.V)  that  it  is 
possible  for  the  pulse  to  overcome  the  band-gap  reflection  and  produce  SIT  in  a  near-resonant 
medium  embedded  in  a  Bragg  reflector.  The  predicted  SIT  propagation  is  a  principally  new 
type  of  a  gap  soliton,  which  does  not  obey  any  of  the  familiar  soliton  equations,  such  as  the 
non-linear  Schrodinger  equation  (NLSE)  or  the  sine-Gordon  equation.  Its  spatio-temporal 
form  and  intensity  dependence  are  shown  here  to  be  distinct  from  the  extensively  —  studied 
gap  solitons  in  Kerr-non-linear  Bragg  reflectors,  which  are  described  by  the  NLSE. 


B.  Results 


In  treatments  of  bidirectional  field  propagation  in  media  with  arbitrary  spatial  distribu¬ 
tion  of  near-resonant  atoms  [31,32],  the  Bloch  equations  for  the  population  inversion  and 
polarization  are  entangled  in  a  fashion  which  leads  to  an  infinite  hierarchy  of  equations  for 
successive  spatial  harmonics.  Here  we  avoid  this  complication  by  confining  the  near-resonant 
two-level  systems  (TLS)  to  layers  much  thinner  than  the  resonant  wavelength,  with  the  same 
periodicity  as  the  dielectric  structure. 

Our  main  idea  has  been  to  try  the  following  phase-modulated  27r— soliton  SIT  solution 
for  the  envelope  of  the  forward  (F)  and  backward  (B)  field 


_  /i  /  1\  exp  [ijanoz/cT,-  At)] 

2/irc  V  u)  °  cosh.  [13 {z / TcCU  —  t)] 


(10) 


where  p,  is  the  transition  dipole  moment,  Tc  is  the  cooperative  (resonant)  absorption  time, 
Aq  is  the  amplitude  of  the  solitary  pulse,  u  is  the  velocity  (normalized  to  c),  no  is  the  mean 
refractive  index  and  A  is  the  field  detuning  from  the  gap  center. 

We  focus  here  on  the  most  illustrative  case,  when  the  TLS  resonance  is  exactly  at  the 
center  of  the  optical  gap.  Then  the  phase  modulation  a,  the  pulse  inverse-width  j3  =  Ao/2 
and  the  detuning  A  are  analytically  obtainable  as  a  function  of  the  group  velocity  cu.  We 
find  that  the  condition  for  SIT  is  that  the  cooperative  absorption  length  cTc/uo  should  be 
shorter  than  the  reflection  (attenuation)  length  at  the  gap  1/k,  i.e.,  that  the  incident  light 
should  be  absorbed  by  the  TLS  before  it  is  reflected  by  the  Bragg  structure.  SIT  is  found  to 
exist  only  on  one  side  of  the  band-gap  center,  depending  on  whether  the  TLS  are  embedded 
in  the  region  of  higher  or  lower  linear  refractive  index  in  the  Bragg  structure.  This  result 
may  be  understood  as  the  addition  of  a  near-resonant  non-linear  refractive  index  to  the 
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modulated  index  of  refraction  of  the  Bragg  structure.  When  this  addition  compensates  the 
linear  modulation,  then  there  is  no  band  gap  and  soliton  propagation  is  possible.  The  soliton 
amplitude  dependence  on  frequency  detuning  from  the  gap  center  (which  coincides  with 
the  TLS  resonance)  is  shown  The  parameters  obtained  from  our  analytical  solutions  fully 
agree  with  those  which  yield  both  forward  and  backward  soliton-like  pulses  in  a  numerical 
simulation  of  Maxwell-Bloch  equations. 

An  adequate  system  for  experimental  observation  of  this  effect  appears  to  be  a  periodic 
array  of  12-nm-thick  GaAs  quantum  wells  (A  =  806nm)  separated  by  A/2  non-resonant 
Al-Ga.As  layers.  Area  density  concentration  cr  ~  10*— 10®  cm"^  of  the  quantum-well  excitons 
yields  ~  lO-^^-lO-^'^s.  A  solitary  pulse  of  <  Ips,  i.e.,  much  shorter  then  the  dephasing 
time  T2  ~  lOps  (at  2°K)  in  this  structure  requires  band-gap  reflection  length  l/«  >  100  A. 

The  salient  advantage  of  the  predicted  near-resonant  gap  soliton  is  stability  with  respect 
to  absorption.  By  contrast,  strong  absorption  is  a  severe  problem  associated  with  a  large 
Kerr  coefficient  required  for  NLSE  gap  solitons  (Appendix  D). 


IV. 2  Theory  of  Quantum  Gap  Solitons 
A.  Background 

By  nature  any  dielectric  material  is  nonlinear.  The  nonlinearity  changes  the  PBG  prop¬ 
erties  of  periodic  structures  to  some  degree.  In  one- dimensional  (ID)  photonic  crystals, 
even  weak  nonlinearity  acts  as  a  large  perturbation  on  these  properties.  In  recent  years  it 
has  been  shown  that  ID  Kerr-nonlinear  photonic  crystals  admit  band-gap  solitary  waves, 
which  are  called  gap  solitons.  Current  studies  correspond  to  classical  gap  solitons.  Since 
quantum  many-body  properties  of  light  can  be  important  in  nonlinear  optical  processes,  we 
have  developed  the  quantum  theory  of  gap  solitons  in  a  ID  Kerr-nonlinear  photonic  crystal 
(Item  5  -  Sec.V'^). 


B.  Results 


The  basic  idea  is  that  incident  photons  with  frequencies  in  a  band  gap  are  scattered 
by  the  nonlinearity  into  the  conduction  and  valence  bands  of  the  photonic  crystal.  The 
effective  Hamiltonian  of  quantum  gap  solitons  is  derived  in  the  two-band  effective-mass 
approximation.  The  eigenstates  of  the  Hamiltonian  are  constructed  exactly  by  Bethe’s 
ansatz  method.  We  find  that  in  a  certain  band  gap  of  the  photonic  crystal  quantum  gap 
solitons  can  be  in  bound  states,  consisting  of  one  or  more  photon  pairs  from  the  valence  and 
conduction  bands.  Such  bound-state  quantum  gap  solitons  are  optical  analogs  of  exciton 
molecules.  Here  we  give  the  2fV-body  bound-state  wave  function  in  the  unnormalized  form 


iV 


.m-Zj-  —  m+Zj^ 


^(^i±,...,^,v±)  cxexp  'ViK- 

m_  - 


/  ^  777-  _  777.  .  I 

X  exp(  -  2^  - - WI Id+  -  D-l 

^  (m_  -  m^)h 


(11) 
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"  j<i 

The  exponential  factors  in  the  bound  state  \^2N  >  can  be  revealed  by  the  dependence  of 
the  intensity-intensity  correlation  function  on  the  separation  rj  of  two  photon  counters 
detecting  the  field  in  the  structure; 

-  j  <  ^2n\^'^{z)S~{z)S^{z +  rj)S~{z -\-'q)\^2N  >  dz  (12) 

where  the  operator  S~{z)  =  4>^{z)  -f  <f>-{z)  is  the  position-dependent  negative- frequency 
field  envelope.  The  2iV-body  bound-state  energy  eigenvalues  are  found  to  have  the  following 
form, 

E2H  =  -Loj2n)  (13) 

Oft  |m_ 

2|m_— m+|  2|m_  — m+l^^  6^^ 

where  -i-  w_)/2  is  the  center  of  the  band  gap  and  the  upper  (lower)  sign  is  chosen 

according  to  that  of  the  Kerr  nonlinearity.  As  seen  from  these  equations,  bound  states  are 
associated  with  discrete  transmission  lines  at  fJa.v  =  cjq  ±  <^2iv  in  the  band  gap. 

As  a  realistic  example,  we  consider  a  periodic  structure  of  alternating  layers  of  GaAs 
crystal  and  linear  dielectic.  The  layer  thicknesses  are  ai(GaAs)  =  0.2338  /im  and  02  =  0.3044 
Hm.  The  refractive  index  of  GaAs  crystal  is  rii  =  3.60  and  the  refractive  index  n2  of  linear 
dielectric  varies  from  1.00  to  3.574.  GaAs  crystal  is  a  self-defocusing  medium  and  has  a 
high  nonlinear  susceptibility  xf  ^  =  —2.5482  x  10“^°  (cm/V)^  for  light  frequencies  below  the 
bandgap  Eg  =  2.1573  x  10^^  s~^.  The  calculation  reveals  that  in  the  second  lowest  band 
gap  of  this  structure  quantum  gap  solitons  can  be  in  bound  states.  The  numerical  results 
reproduce  the  multi-exponential  fall-off  of  the  intensity-intensity  correlation  function 
with  the  detector  separation  77  and  the  variation  of  the  optical  exciton  frequency  Q,2N  with 
the  relative  refractive  index  ni/n2,  for  different  quantum  numbers  N. 

In  conclusion,  bound-state  quantum  gap  solitons  will  manifest  themselves  via  distinct 
dependence  of  the  intensity-intensity  correlation  function  on  the  detector  separation  and  via 
transmission  resonances  at  band-gap  frequencies  (Appendix  E). 


IV. 3  Superluminal  Delays  of  Coherent  Pulses  in  PBG  structures 

A.  Introduction 

As  shown  by  a  two-photon  interference  experiment  [33],  a  photon  that  has  tunneled  as 
an  evanescent  wavepacket  through  a  PBG  structure  (dielectric- mirror)  appears  to  have  been 
delayed  significantly  less  than  its  ’’twin”  photon  that  has  traversed  the  same  distance  in 
vacuum.  Such  a  delay  has  been  interpreted  as  signifying  ’’superluminal”  (faster-than-light) 
barrier- traversal  time.  Similar  ’’superluminal”  time  delay  in  tunneling  through  a  dielectric 
mirror  has  now  been  measured  in  a  classical  two-pulse  interference  experiment  [34].  The 
latter  experiment  has  also  revealed  a  remarkable  feature,  namely,  that  the  temporal  width 
of  the  transmitted  wavepacket  is  strongly  narrowed  down. 
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It  is  always  possible  to  trace  numerically  the  evanescent  wavep^cket  evolution  and  com¬ 
pare  its  features  with  different  definitions  of  barrier  traversal  times  [35-40]  (see  below). 
Nevertheless,  the  mechanism  of  superluminal  time  delays  is  still  obscure  [33]  and  regarded 
as  a  ’’poorly  resolved  mystery”  [40].  A  commonly  invoked  notion  is  that  this  mechanism  is 
spectral  reshaping  (filtering)  of  the  transmitted  wavepacket  by  dispersion.  Yet  why  should 
spectral  reshaping  necessarily  yield  superluminal  delays  of  EM  pulses  in  non-absorbing  struc¬ 
tures,  after  propagating  in  (dispersionless)  vacuum?  Is  there  a  common  mechanism  for 
superluminal  time  delays  and  wavepacket  narrowing,  which  applies  to  both  EM  pulses  in 
dielectric  structures  and  relativistic  massive  particles  in  potential  barriers  [40]?  How  is 
causality  compatible  with  superluminal  transmission,  particularly  in  the  single-photon  case 
[33]? 

We  purport  to  show  (Item  3  "■  Sec.V)  that  the  above  questions  can  only  be  answered  by 
a  universal  description  of  the  temporal  wavepacket  transmission  as  interference  between  its 
causally-propagating  consecutive  components.  Our  description  reveals,  for  the  first  time,  the 
key  role  of  phase  coherence  in  tunneling,  by  demonstrating  its  dependence  on  the  coherence 
time  (phase  randomization)  of  the  wavepacket.  An  important  corrolary  is  that  superluminal 
time  delays  can  occur  also  in  allowed  propagation,  namely,  propagation  which  can  only  be 
described  by  real  wavevectors  (e.g.,  in  Fabry-Perot  structures)  and  not  only  in  evanescent- 
wave  tunneling,  where  complex  wavevectors  can  be  employed  (e.g.,  in  photonic  band-gap 
structures). 


B.  Results 

Our  theory  has  demonstrated,  for  the  first  time,  that  the  universal  mechanism  of  predomi¬ 
nantly  destructive  interference  between  accessible  causal  paths  is  responsible  for  transmission 
attenuation,  superluminal  delay  times  and  wavepacket  narrowing.  Two  other  characteristics 
of  evanescent  waves  ,namely,  exponential  attenuation  and  traversal-length  independence  of 
the  mean  traversal  time  of  the  structure  length  [34]  can  also  be  explained  in  terms  of  this 
universal  mechanism.  This  theory  overcomes  the  limitations  of  previous  approaches,  since 
it  applies  to  arbitrary  pulse  shapes,  widths  and  coherence  times,  and  explicitly  reveals  the 
causal  nature  of  their  transmission.  The  understanding  provided  by  this  theory  may  open 
new  perspectives  in  the  design  of  the  velocity,  intensity  and  shape  of  transmitted  pulses,  by 
manipulating  the  phase  delays  along  the  accessible  paths  in  the  medium  (Appendix  F). 
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Appendix  A: 

Non-Classical  Field  States  via  Periodically  Non- Adiabatic 

transitions 

B.  Sherman,  G.  Kurizki  and  A.  Kozhekin 
V.  Akulin^  and  A.  Levine^ 

Chemical  Physics  Department,  Weizmann  Institute  of  Science,  Rehovot  76100  Israel 

A  novel  mechanism  for  generating  nonclassical  field  states  has  been  dis¬ 
covered,  based  on  non-adiabatic  periodic  modulation  of  the  parameters  gov¬ 
erning  the  field-atom  coupling. 


The  resonant  interaction  of  a  quantized  field  with  a  two-level  atom  has  been  intensively 
studied  in  single-mode  resonators  where  the  atomic  motion  plays  no  role,  since  the  field 
is  taken  to  be  spatially  uniform  in  the  beam  direction.  The  field-atom  interaction  then 
conforms  to  the  fundamental  Jaynes-Cummings  model  [1,2]  (JCM). 

Our  purpose  here  is  to  show  that  atomic  motion  along  a  cavity  mode  can  give  rise  to 
new  QED  effects  that  originate  from  the  spatial  character  of  the  field.  These  effects  can  be 
revealed  by  propagating  an  atomic  beam  through  a  cavity  along  its  axis  or  through  a  defect 
in  a  photonic  crystal  [3]. 

It  has  been  shown  that  periodic  temporal  modulation  of  the  resonant  field- atom  coupling 
for  an  atom  moving  along  a  standing-wave  mode  can  drastically  modify  the  prominent 
dynamical  features  of  the  JCM  [1,2],  such  as  (a)  the  evolution  of  ’’Schrodinger  cats”  and  their 
photon  statistics  [4];  (b)  the  corresponding  oscillations  of  the  atomic  population  inversion 
[5].  Hear  we  show  that  near-resonant  oscillatory  modulation  of  the  coupling  in  the  same 
system  gives  rise  to  even  more  dramatic  modifications  of  atomic  and  field  dynamics.  In 
addition  to  their  inherent  novelty  and  interest,  these  modified  dynamical  features  are  shown 
to  allow  much  better  control  than  their  JCM  counterparts  over  the  generation  of  nonclassical 
field  states. 

The  system  under  consideration  involves  the  classical  translational  motion  of  the  atom, 
the  internal  atomic  states,  and  the  quantized  EM  field.  The  atomic  internal  states  become 
strongly  entangled  with  the  quantized  field  state  in  the  course  of  their  interaction.  The 
resulting  combined  field- atom  system  has  eigenenergies  and  eigenstates  (’’dressed  states”). 
The  dressed-state  energies  depend  on  the  field  amplitude  and  therefore  on  the  atomic  lo¬ 
cation  within  the  field  mode.  In  the  limit  of  slow  atomic  motion  along  a  standing- wave 
field  mode,  the  field-atom  system  adiabatically  follows  the  dressed  state  associated  with  the 
instantaneous  atomic  position.  However,  the  system  can  still  undergo  non-adiabatic  transi¬ 
tions  from  one  dressed  state  to  another,  with  exponentially  small  probability.  In  this  case, 
the  kinetic  energy  of  the  classically-moving  atom  is  transferred  to  the  internal  quantum 
degrees  of  freedom  of  the  combined  field-atom  system.  Such  a  non-adiabatic  transition  is 
very  unlikely  for  sufficiently  slow  atoms,  and  therefore  it  occurs  near  the  field  nodes,  where 
the  energy  required  for  such  a  transition  is  minimal  (see  Fig.l).  Since  the  transition  can 
occur  at  any  of  the  nodes  along  the  atomic  path,  the  amplitudes  of  the  transitions  at  dif¬ 
ferent  nodes  interfere.  The  phases  of  these  interfering  amplitudes  are  very  sensitive  to  the 
number  of  photons  in  the  quantized  field.  Hence,  only  for  several  photon-number  states  the 
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amplitudes  may  interfere  constructively^  significantly  contributing  to  the’ probability  of  the 
non-adiabatic  transition.  This  implies  that  non-adiabatic  transitions  are  strongly  correlated 
with  certain  photon-number  states.  Therefore,  a  measurement  that  finds  an  atom-field  sys¬ 
tem,  which  was  initially  prepared  in  the  lower  dressed  state,  in  the  upper  dressed  state  (which 
almost  coincides  with  the  atomic  ground  state  in  the  case  of  large  detuning  A  =  cua  —  cj), 
projects  the  initial  field  distribution  on  these  number  states  out.  The  final  field  statistics 
is  the  convolution  of  the  initial  one  with  the  weights  of  different  photon  numbers  in  the 
non-adiabatic  transition  probability. 

The  Hamiltonian  for  a  two-level  atom  with  dipole-transition  frequency  moving 
through  the  electromagnetic  field  of  a  spatially  periodic  mode  with  frequency  u  reads 

H  =  -f-  hua^ ah^[f{t)\{a +  d^d_).  (1) 

Here  d+((7_)  and  a^(a)  are  the  atomic  and  field  raising  (lowering)  operators,  respectively, 
whereas  the  coupling  strength  fi(r)  depends  on  the  amplitude  of  the  field  mode  at  a  given 
point 


hQ,[f{t)]  =  (^w/2)^^^  (/?  •  ex)£  {z  =  vt  +  fx)  (2) 

.  where  p  is  the  atomic  dipole  moment  and  t\  is  the  mode  polarization. 

We  consider  classical  translational  motion  of  an  atom  with  the  constant  velocity  v  along 
the  nodes  of  the  mode  in  a  structure  of  length  1.  We  assume  that  the  kinetic  energy  of  the 
atom  is  much  bigger  than  the  energy  of  the  non-adiabatic  transition,  such  that  the  dispersion 
of  the  wave  packet  due  to  the  velocity  difference  after  and  before  the  the  transition  are 
much  smaller  than  the  atom's  size,  that  is  hA./mv^  ^  Ratll-  This  condition  is  satisfied, 
for  example,  for  a  cesium  atom  with  velocities  of  10^ cm /s  and  non-adiabatic  transition 
frequency  A  of  moving  through  a  resonator  with  a  length  /  of  1cm. 

For  an  atom  moving  in  a  sinusoidally  periodic  mode  with  period  A  the  amplitude  S{z) 
can  be  chosen  to  be  real,  here  we  can  write 

f2[2:(f)]  =  Qo  sin(27rf/T),  (3) 


with  T  —  X/v. 

The  Hamiltonian  (1)  has  two  adiabatic  energy  terms 

e±(f)  =  hoj{n  +  1)  ±  ^v^A2(t)/4  +  n^{t){n  +  1).  (4) 

In  the  limit  of  vanishingly  small  atomic  velocity,  the  evolution  of  e^{t)  is  adiabatic,  i.e,,  if  ini¬ 
tially  the  atom  was  in  the  eigenstate  with  energy  e^(0),  then,  after  time  t,  the  eigenstate  ac¬ 
quires  an  adiabatic  phase  factor  exp[— i  ef{£)dt'/h]  Thermal  velocity  v  of  the  off-resonant 
atom  results  in  weak  non-adiabaticity,  if  the  Doppler  frequency  uj^  =  27r/T  =  (27r/A)u,  is 
small  in  comparison  with  the  atom-field  detuning  A  =  —  w,  A  >  LUp.  Non-adiabatic 

transitions  are  possible  only  in  the  vicinity  of  the  points  of  the  energies  closest  approach 
[6],  namely,  where  £l[t)  =  0  (see  Fig.l).  In  order  to  find  the  probability  amplitude  of  such 
a  transition,  one  can  in  principle  use  the  first  order  approximation  in  the  non-adiabatic 
coupling  [7].  However,  since  the  result  is  exponentially  small,  the  first  order  approximation 
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in  the  non-adiabatic  coupling  yields  incorrect  result  [8].  A  more  precise  expression  for  the 
amplitude  of  the  transition  Un  is  derived  by  Dykhne  [8]. 


On  =  exp 


Im  dtyj /S?  +  4fi^(i)(n  +  1)  , 


(5) 


where  tc  is  the  branching  point  of  the  integrand  (the  points  of  levels  crossing)  in  the  complex 
time  plane,  that  is  closest  to  the  Re  t  axis. 

In  this  approximation  the  non-adiabatic  transitions  are  sudden  jumps  at  the  nodes  of  the 
standing  wave.  The  phase  difference  for  transitions  at  two  adjacent  transition  space-points 
reads 

dt\J -f-  4flo(^  +  1)  sin^tufj^,  (6) 

and  results  in  the  following  net  amplitude  of  the  non-adiabatic  transition  after  K  nodes  the 
atom  has  passed 


K 

E 

;=o 


il(pn 


sin  Kipn 

—  i 

sm  ifn 


(7) 


Since  the  probability  of  the  non-adiabatic  transition  is  small,  we  neglect  processes  with  more 
than  one  net  jump  (the  next  nonvanishing  order  of  approximation  would  involve  triple  jumps 
—  from  the  lower  to  the  upper  state,  then  down  to  the  lower  state  and  finally  back  to  the 
upper  state). 

We  now  consider  the  field  evolution  and  assume  that  the  atom  has  been  initially  prepared 
in  its  ground  state  l^'),  that  is 

|■^(0))  =  Z!cn|n-b  l)lg).  (8) 


The  atom  is  taken  to  start  and  finish  its  motion  at  the  nodes,  where  =  0  and  the  ground 
and  excited  states  are  the  lower  and  upper  dressed  states,  respectively.  Hence,  if  after  the 
interaction  with  the  field  the  atom  is  found  in  the  upper  state,  this  means  that  the  non- 
adiabatic  transition  has  definitely  happened.  The  field-atom  final  wavefunction,  reduced  by 
the  conditional  measurement  of  atomic  excitation  is 


\^{t  =  KT))  =  NJ2  CnAne'^‘^^^\n)\e),  (9) 

n 

where  N  is  the  normalization  factor.  We  see  from  (9)  that  the  resulting  field  statistics  is  the 
convolution  of  the  transmission  function  and  the  initial  photon  distribution  c„.  This  can 
result  in  selecting  ("filtering”)  out  either  a  single  or  several  number  state  components,  with 
the  corresponding  phases,  thus  generating  a  strongly  non-classical  state.  We  can  adjust 
over  a  wide  range  of  possibilities  via  the  following  dimensionless  parameters:  p  =  A^/dffo, 
^  =  fio/wi?  and  the  number  of  nodes  of  the  standing  wave  K.  The  net  amplitude  A„  has 
narrow  peaks  corresponding  to  the  constructive  interference  condition  on  the  adiabatic  phase 
difference  (6). 

(/?„  =  7rm,  m  =  1,2,3....  _  (10) 
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Ijnder  this  condition,  the  multichannel  constructive  interference 'enhances  the  transition 
probability  by  K^.  If  only  one  photon  number  produces  the  peak,  which  contributes  to  the 
convolution  (9)  with  the  initial  field  distribution,  then  the  final  field  state  becomes  nearly  a 
Fock  state  (see  Fig. 2)  after  the  conditional  measurement  of  the  atomic  excitation.  Thus  a 
presecelected  Fock-state  is  generated. 

Consider  now  the  opposite  case  of  wide  range  of  n  satisfying  the  constructive  interference 
condition  (10)  and,  therefore,  contributing  to  the  convolution  (9).  We  assume  that;  (i)  the 
adiabatic  phase  difference  ipn  ~  (ii)  a„  do  not  vary  strongly  with  n.  This  means  that 
one  gets  the  convolution  of  the  initial  filed  statistics  with  the  sine  of  -\/n,  which  is  typical  to 
one-photon  JCM  [1,2].  As  it  was  shown  [9-11]  for  a  field,  initially  prepared  in  a  quasiclassical 
coherent  state  a  conditional  measurement  leads  to  the  formation  of  a  macroscopic  quantum 
superposition  of  two  nearly- coherent  states  with  identical  mean  amplitudes  (same  as  the  ini¬ 
tial  one)  and  different  mean  phases  (’’Schrodinger  cat”).  The  corresponding  Q-function  and 
the  phase  distribution  F'(d),  described  in  the  Pegg-Barnett  formalism  [12]  by  the  function 
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(11) 


are  shown  in  Fig. 3. 

The  intermediate  regime,  when  several  number  states  are  singled  out  from  the  initial 
distribution  by  the  constructive  interference  and  the  conditional  measurement  (CM)  of  the 
atom  in  |e),  results  in  the  generation  of  amplitude  squeezed  states  similar  to  the  ones  predicted 
by  \amamoto  et  al.  for  strong  self-squeezing  in  Kerr  media  [13]  (  see  Fig. 4). 

The  weakly-nonadiabatic  regime  considered  above  imposes  two  major  limitations  on 
non-classical  field-state  preparation  via  constructive  interference  of  nonadiabatic  transitions 
followed  by  a  CM  of  the  atom:  (i)  the  CM  success  probability  is  typically  low  (  1%);  (ii)  the 

resulting  Fock-state  amplitudes  are  controlled  by  very  few  parameters,  which  restricts  the 
type  of  field  state  that  can  be  prepared.  In  order  to  overcome  these  limitations,  we  adopt 
the  model  of  temporal  Kronig-Penney  modulation,  e.i.  ,  periodically  alternating  abrupt 
jumps  of  the  field  detuning  A  or  coupling  constant  (Fig. 5).  Then  the  nonadiabatic  transition 
probability  is  large  in  every  cycle  of  the  modulation.  This  type  of  modulation  can  be  achieved 
in  photonic  band  structures  [3],  when  the  mode  frequency  is  near  the  forbidden  gap  edge. 
Under  the  condition  of  constant  atomic  velocity,  the  atom-field  coupling  is  modulated  with 
period  T  =  tp  +  t^),  where  tf  is  the  duration  of  the  zero-coupling  stage  (A  ^  c<2£),flo)j  and 
to  is  the  duration  of  the  ‘dressed’  evolution  (A  ~  f^o)-  The  2x2  evolution  matrix  of  the 
first  stage  is  diagonal,  it  describes  the  free  oscillations  of  atom+field  states  with  appropriate 
energies 


Ufitp)  =  exp[— i(u;d^d  -f  (w  —  /S.)crzl2)tF]  =  *^a^F/2) 


exp{iua.tF  12)  J 

(12) 


The  system  evolution  at  the  ‘dressed’  part  of  the  cycle  is  described  by 

UoiiD)  =  exp{— i[cju^d  -f  (cu  —  A)<J3/2  +  flo(d-+d  -f-  a^a^)]to}-  (13) 

The  later  expression  can  be  written  in  the  matrix  form  as  follows 
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Ud{^d)  — 


'  g-tu/a<D/2(cos 'yi£)  —  ^sin^ 

atD/2 


_ji2a4ig‘“at 
7 


(14) 


Here  /3  =  +  AV4,  7  =  ^^^(ata  +  1)  +  A2/4 

The  abrupt  transition  between  the  two  stages  allows  to  use  the  sudden  approximation, 
i.e.,  to  assume  that  the  system’s  state  remains  unchanged  during  the  switching.  Thus  the 
single-cycle  evolution  operator  is 

U{T)  =  UF{iF)UD{tD).  (15) 


The  entangled  field-atom  state  after  K  cycles  of  modulation  is  given  by 


m  =  KT))  =  U'^(T)\i,m  =  [UF(tr)fMtD)f\i>m. 


(16) 


The  field  state  projected  by  successful  measurement  of  the  atom  in  |e)  after  K  cycles 
has  the  form,  analogous  to  expressions  (7-9),  with  ipn- 


(fn  =  arccos  <  cos  {/^nto)  cos  {AtF/2) 


2^n 


cos  {^nto)  cos  (Atir/2) 


(17) 


It  is  possible  to  design  more  complicated  modulation  schemes,  in  which  each  cycle  con¬ 
sist  of  several  unequal  intervals  with  different  detunings.  We  can  thereby  control  enough 
parameters  to  prepare  a  superposition  of  several  Fock  states  with  arbitrary  amplitudes  from 
an  initial  coherent  state,  by  the  same  as  described  above.  In  Fig. 6  we  show  that  (a)  superpo¬ 
sition  of  3  Fock  states  with  arbitrary  pre- determined  amplitudes  can  be  obtained  from  initial 
coherent  states,  using  modulation  with  3  different  detunings  and  time  intervals  within  each 
cycle;  (b)  The  CM  success  probability  is  much  larger  than  the  squared  project  of  the  initial 
onto  the  final  state. 
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Appendix  B: 

Lasing  Without  Inversion  via  Spontaneously  Formed  Coherence 

in  Photonic  Band  Gaps 

A.G.  Kofman  and  G.  Kurizki 

Chemical  Physics  Department,  Weizmarm  Institute  of  Science,  Rehovot  76100,  Israel 
Assume  that  the  transition  between  a  ground  I5)  and  an  excited  |e)  levels  of  the  atom  is 
near  resonant  to  a  photonic  band.  The  strength  of  the  atom-band  coupling  is  characterized 
by  the  function 


^(‘^)  =  Z!  (1) 

A 

where  k\  is  the  matrix  element  of  the  coupling  between  the  atom  and  the  mode  A  with 
the  frequency  cja-  The  band  is  assumed  to  be  well  separated  from  other  bands,  so  that  the 
influence  of  the  latter  can  be  neglected.  Moreover,  the  band  is  assumed  to  be  sufficiently 
dense  and  narrow,  i.e.. 


Gra  >  Aj,  |A„,|, 


(2) 


where  Gm  and  A;,  are  the  maximum  value  and  the  width  of  the  function  G{(jj)  respectively. 
Am  =  LOa  —  <^Tn,  is  the  Center  of  gravity  of  the  band, 

Um  =  ~r  duJuG{u}),  (3) 

JuJxj 

(jju  and  cjl  are  the  lower  and  upper  band  edges  respectively,  and 

r^L 

duG{u).  (4) 

JlJu 

Then  the  level  e  can  be  shown  [1]  to  split  to  two  levels  1  and  2,  described  by  the  normalized 
wavefunctions. 


\i’j)  = 


®)0a)  +  Z 


Uj  -  ujx 


-I9.U) 


0  =  1,2). 


(5) 


The  energies  hui^2  and  the  amplitudes  of  the  states  1  and  2  can  be  expressed  through 
the  moments  of  G{u)).  In  particular,  the  level  separation  to  first  approximation  is 

U2  —  LOi  =  -|-  A^  =  (6) 

where  Q/j  is  the  vacuum  Rabi  frequency.  An  excited  atom  injected  at  t  =  0  into  the  PBS 
in  the  case  (2)  undergoes  an  insignificantly  small  spontaneous  decay  over  a  time  of  order 
l/rifi,  after  which  the  excited-state  wavefunction  becomes  a  coherent  superposition  of  the 
states  IV’i)  and  |■02)  [1], 

l^(t))  =  -b  -  (7) 
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If  there  is  a  nonzero  matrix  element  of  the  dipole  moment,  dse,  for  the  transition  between  the 
level  e  and  another  level  3,  it  then  follows  from  (5)  that  the  dipole  moment  matrix  elements 
for  the  transitions  1-3  and  2-3  are  given  by 


^31  —  V^dse,  dz2  =  x/c^dse.  (8) 

For  the  narrow-band,  near-resonant  case,  iA,  which  is  of  special  interest  here,  one 

gets  in  the  first  approximation 


1  rr 

Cl  =  C2  =  CUi,2  =  - - -  T  V^.  (9) 

The  density  matrix  elements  of  the  system  consisiting  of  the  levels  1  and  2  are  [cf.  Eqs.  (7) 
and  (9)] 


piiit)  ^  p22{i)  =  2’  ^12(0  =  p2i{t)  ^  ^ — •  (10) 

The  coherence  between  the  levels  1  and  2  can  be  utilized  to  produce  lasing  without 
inversion  [2].  Following  Refs.  [3-5],  consider  a  A-system.  The  field 

E(t)  =  Eqc  +  c.c.  (11) 

couples  the  upper  level  3  with  the  lower  levels  1  and  2,  The  levels  1,  2,  and  3  have  have  the 
reciprocal  lifetimes  respectively  71,  72,  and  73  =  731+732  +  73:  where  7^^,  7^2:  ^nd  7'  are  the 
spontaneous  rates  of  transitions  from  the  level  3  to  the  levels  1,  2,  and  all  others  respectively. 
The  levels  1  and  2  can  have  finite  lifetimes  due  to  one  or  several  of  the  following  factors: 
spontaneous  transitions  to  levels  other  than  the  ground  state,  finite  time  of  presence  of  the 
atom  in  the  structure,  ’’background"’  of  modes  (outside  the  narrow  band  discussed  above, 
which  causes  the  splitting  of  level  e),  and  collision-induced  transitions.  The  open  3-level 
system  consisting  of  the  levels  1,  2,  and  3  is  described  here  by  an  unnormalized  density 
matrix  i?(t),  whose  diagonal  elements  yield  the  number  of  atoms  in  the  laser  cavity  in  the 
respective  level.  In  the  rotating  wave  approximation,  the  rotating- frame  matrix  elements  of 
R{t)  obey  the  equations, 

^11  =  —  V^lRlz)  —  JlRu  +  731^33  +  I+ii, 

R22  —  ~^(^3-^32  ^2-^23)  ““  'I2R22  +  732^33  +  1+225 

^33  —  —{{V^iRis  +  VS2R23  ~  VlsRsi  —  V23R32)  ~  I3R33  +  1+335 
i?12  ==  {i(^2l  j2l)Rl2  +  ^(^2^13  ~  ^13-^32)  +  1+125 

Ri3  =  (2^31  —  73i)i?i3  +  ^V23^12  ~  ‘^1^3(^33  ~  ^11)5 

R23  =  (^*^32  —  732)^23  +  ^+13^21  “  ^^3(^33  ~~  -^22)*  (12) 

Here  Asj  =  cj,  Vsj  z=:  =  —d^jEojU  [j  =  1,2),  is  the  decay  rate  of  the  nondiagonal 

matrix  element  VF33  is  the  pumping  rate  of  atoms  to  level  3,  Wij  =  gp?-  (z,j  ~  1,2),  q 
being  the  rate  at  which  atoms  form  a  coherent  superoisition  of  the  states  1  and  2  (by 
spontaneous  decay  of  |e)),  described  by  the  normalized  density  matrix  p?-. 

The  complex  susceptibility  x  is  expressed  through  the  solution  of  Eqs.  (12)  by 


22 


(13) 


X  -  X^  -  ^x"  —  y^[.dz\Rlz  +  dz2R23) 


(in  the  Gaussian  system  of  units),  where  V  is  the  structure  volume.  The  steady-state  solution 
of  Eqs.  (12)  yields  in  the  first  order  in  the  field, 

_  f  \diz\^ Niz  +  f  \d23\'^N23+J^\ 

^  "  nv[  'y3i-  ^ A31  732  -  ^■A32  ;  ’  ^  ’ 

where  Niz  =  -  i?33  {i  =  1,2),  is  the  steady-state  value  of  Rij  in  the  absence  of  the 

field. 


^(0)  _  )T33  ^(0)  ^ 


VE„  +  -^W33  (i  =  l,2), 

73  / 


e  j  j  d(0)  dzid23^12 

f  -  dzid23Ri2  =  - : - • 

721  -  2^21 


Eqs.  (13)  and  (14)  yield 


Re  (  ^  _|_  l^23p-^^23  +  / 

V  731  ~  'iAsi  732  +  f  A32 


The  necessary  condition  for  lasing  is  x"  <  0-  Below  we  consider  some  situations  where  lasing 
is  possible  in  the  absence  of  inversion,  i.e.,  for  Niz  >  0  and  Nzz  >  0. 

Let 

731  =  732  (18) 


A31  =  —  A32,  i.e.,  ui  —  (wsi  -l-a;32)/2. 


Then  for 


arg  f  =  IT  —  <j),  where  tan  <j)  =  7; - ,  (20) 

2731 

the  terms  proportional  to  /  in  Eq.  (17)  are  negative,  yielding  with  the  account  of  (16)  and 
(17) 

/,  _  47r  2|d3xd32VEi2|  _  73l(|d31piVl3  +  |d32|^iV23)  ^2]^^ 

xlllx  +  ^21/ 4  .  \j 721  +  W21  \lll\  +  ‘*^21/4 


Consider  first  the  case  of  a  very  small  separation  of  the  levels  1  and  2,  0221  721,  73i- 

In  this  case  <l>  ~  !jJ2\/{‘^13\)  1  in  Eq.  (20).  Then 


23 


(22) 


x"  = 


47r 


^^731 

+ 


33 


1  -  |<i3.P  + 

V  7i  y 


(’-s 


\d: 


32 


Msi  p^ll  ^  |(i32pVV22  ^\d3ld32Wi2\ 


7l  72  721 

The  second  bracketed  term  here  can  be  shown,  with  the  help  of  inequalities  721  >  (71  +72)/2 
and  /}°2  ^  \/puP22j  be  nonnegative.  However,  it  can  be  made  very  small.  In  particular, 


it  vanishes  if 


lAi 


Pii 


-0 

P221 


1^31 1  —  1^32], 


7i2  =  7i 


72- 


Then 


X 


8;r|d; 


31 


hv-f3i 


1 


731  ~b  732 

271  , 


(23) 

(24) 

(25) 


If  731  +  732  <  271  the  lasing  is  possible  for  very  small  number  of  atoms  at  the  levels  1  and 
2,  irrespective  of  the  populations  of  the  levels  1  and  2. 

Consider  now  the  opposite  case  of  a  large  separation  of  the  levels  1  and  2,  u>2i  '>  721, 2731. 
In  this  case  «  ir/2  —  273i/c<;2i  ~  7r/2  in  Eq.  (20).  Assuming  the  validity  of  (23)  which  is 
justified  in' the  above  closely  resonant  case  of  the  level  splitting  in  a  PBS  [cf.  (8)-(10)]  and 
assuming  also  that 


7i  =  72, 


Eq.  (21)  becomes 


// 

X  = 


ISiryildsil 


731  ~  71/2, 

(26) 

.(0)  Pi  ^  d(0) 

1  273.)^'+ 

(27) 

If  731  =  (71  +73)72,  the  second  inequality  in  (26)  is  equivalent  to  73  <C  71  and  (27)  becomes 


X 


*'33 
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ll 


Equations  (27)  and  (28)  show  that  due  to  a  coherent  superposition  of  the  states  1  and  2 
only  a  very  small  fraction  of  the  atoms  in  these  states  absorb  the  field  E{t)  and  therefore 
can  yield  gain  when  the  number  of  the  atoms  at  the  level  3  is  much  less  than  the  number  of 
the  atoms  at  the  lower  levels  1  and  2. 

This  work  has  been  supported  by  USARDSG. 
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Appendix  C: 

Photon  Exchange  by  Atom  Pairs  in  Resonators 
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A  system  of  two  identical  atoms  sharing  a  photon  with  a  high-Q  resonator 
mode  is  studied  by  a  non-perturbative  formalism.  As  opposed  to  previously 
considered  models,  here  the  interatom  coupling  is  shown  to  residt  from  com¬ 
peting  effects  of  vacuum  Rabi  splitting  and  photon  exchange  via  off-resonant 
modes.  Strong  suppression  of  interatom  excitation  transfer  is  predicted  at 
both  near-  and  far-zone  separations. 

PACS  numbers:  42.50.Fx,  42.50.Md 

The  physical  importance  of  a  system  of  two  resonant  atoms  sharing  a  photon  is  evident 
from  the  numerous  studies  of  its  various  aspects;  causality  and  retardation  of  the  interatomic 
photon  exchange  [1];  atomic  level  shifts  due  to  the  photon  exchange,  which  are  identical  with 
the  resonant  dipole-dipole  interaction  (RDDI)  [2,3];  oscillatory  cooperative  decay  in  diatom 
dissociation  [4];  spectra  [5]  and  squeezing  [6]  of  two-atom  resonance  fluorescence.  The  above 
literature  has  dealt  with  atoms  in  open  space,  whereas  studies  of  cooperative  atomic  effects 
in  resonant  cavities  have  been  extensions  of  the  Tavis- Cummings  model  [7].  This  model 
assumes  many  atoms  identically  coupled  to  a  single  mode  and  ignores  the  symmetry- breaking 
dipole-dipole  effects,  which  are  important  at  near-zone  separations  [3,5]. 

This  Rapid  Communication  is  aimed  at  gaining  new  insight  into  systems  which  are  now 
becoming  experimentally  realizable  [8]:  a  pair  of  identical  two-level  atoms  (or  excitons), 
sharing  a  photon  with  one  or  many  high-Q  modes  in  a  resonator.  To  this  end,  we  develop  a 
formalism  capable  of  treating  two  atoms  coupled  to  a  field  with  an  arbitrary  mode-density 
spectrum  in  a  non-perturbative  fashion.  Such  a  formalism  is  necessary  because  the  standard 
perturbative  treatment  of  two-atom  coupling,  to  second  order  in  the  field  [2],  is  inadequate 
for  a  reservoir  whose  mode-density  spectrum  does  not  vary  smoothly,  as  demonstrated  al¬ 
ready  by  our  comprehensive  theory  of  a  single  atom  coupled  to  such  a  reservoir  [9].  We 
predict  drastic  modifications  (as  compared  to  previously  studied  models)  of  the  energy  lev¬ 
els  and  interatom  excitation  transfer  in  such  systems,  both  at  near-zone  and  fax-zone  atomic 
separations.  These  modifications  are  due  to  interfering  (competing)  effects  of  strong  atomic 
coupling  to  high-Q  modes  (vacuum  Rabi  splitting  [7])  and  interatom  photon  exchange  via 
all  other  modes  (RDDI  [2,3]). 

Here  we  concentrate  [10]  on  a  system  of  atoms  A  and  B  that  are  nearly-resonant  with 
a  cavity  mode  or  band  of  degenerate  modes  (as  in  spherical  resonators).  The  effective 
Hamiltonian,  which  can  be  derived  from  first  principles,  has  then  the  following  second- 
quantized  form  in  the  rotating-wave  approximation  (RWA) 

H  =  h{u)A\^A){eA\+<^B\eB){^B\) 


+  hM{\eA){eB\  +  \eB){eA\) +  fi'^ioxaxa\ 

A 

+  ^'^{f^Axax\&A){gA\  +  >^Bxax\eB){gB\  +  H.C.).  (1) 

A 

Here  |e4(B))  and  |5'a(B))  are  the  excited  and  ground  atomic  states.  The  frequency  ujx  and 
annihilation  operator  ax  pertain  to  a  near-resonant  field  mode  labeled  by  A  whose  position- 
dependent  dipolar  coupling  to  atom  A  or  H  is  given  by  ka{b)]  ^'^a{B)  are  the  excited-state 
energies  including  the  spatially-dependent  shifts  due  to  single-atom  interactions  with  the 
field  reservoir  outside  the  near-resonant  narrow  line  (band)  and  M  is  the  matrix  element 
of  the  interatom  interaction  via  the  same  reservoir  calculated  by  second-order  perturbation 
theory  [2,3].  Both  fiiOA(B)  and  M  may  have  imaginary  parts  (see  below). 

The  time-dependent  wavefunction  for  the  system  of  a  single  photon  shared  by  the  field 
and  atoms  can  be  written  as 


I'f6(f))  —  OCAb{t)\eAgB,  {Oa})  +  OiBb{t)\eBgA,  {Oa}) 

+  1a),  (2) 

A 

where  |{0a})  and  |1a)  are  the  vacuum  and  A-mode  single-photon  states,  respectively,  whereas 
^  =  A,  B  denotes  which  atom  is  excited  initially.  The  Schrodinger  equation  yields  the 
following  exact  Laplace-transform  solutions  (within  the  RWA)  for  the  excited-state  S'-matrix 

d(s)  =  i}-i(s)6^(s).  (3) 

Here  q:(s)  and  U{s)  are  2  x  2  matrices  in  the  basis  of  \eA(B)gB{A))-  The  diagonal  elements 
of  U{s)  are  Uaa(bb)  =  s  +  i'-^B(A)  +  *Tb(^)(s),  where 


_  [  Ga(b){^)<Bj 
is  —  ui  ’ 
Ga(B){<^)  -  -  wa). 

A 


JA{B)is)  -  j 


(4) 


Here  and  henceforth  the  integration  is  performed  only  over  the  narrow  band  of  near-resonant 
modes.  The  off-diagonal  elements  are  Uab{ba){s)  =  —iM  —  iJAB{BA){s)i  where 


J  IS  --  u 

Gab(ba)M  =  ^a(b)xk*b(a)xH^  -  ^x)-  (5) 


The  resolvent  (denominator)  of  Eq.  (3),  £>(3),  is  the  determinant  of  the  matrix  B(s). 
The  roots  of  the  equation  D(s)  =  0,  which  correspond  to  the  levels  (eigenvalues)  of  the 
system,  can  differ  strongly  from  the  standard  perturbative  solutions,  as  shown  below.  Our 
subsequent  analysis  rests  on  the  following  assumptions,  (i)  The  transition  frequencies  and 
decay  rates  of  atoms  A  and  B  are  equal,  which  is  true  for  identical  atoms  far  from  cavity 
walls  or  dielectric-layer  interfaces,  (ii)  The  transition  dipole  moments  are  real  and  the 
mode  functions  are  standing  waves,  i.e.,  also  real;  then  Gab(‘^)  =  Gba(<ja)  is  real,  yielding 
Uab{3)  =  Uba{s),  whence  Q!.4b(^)  =  aBA{t)- 
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The  case  of  weak  atom-reservoir  coupling  obtains  if  the  near-resonant  line  (band)  is 
sufficiently  broad  and  not  sharp  (e.g,  a  single  mode  broadened  due  to  its  finite  (^-factor). 
By  suitable  extension  of  the  validity  conditions  for  the  Weisskopf-Wigner  approximation  for 
single-atom  emission  [9],  we  find  that  this  case  holds  if  the  mode-density  spectrum  is  smooth 
enough,  so  that 


dJA{B) 

duiA 


dGA(B){^A) 

(ko 


Ga(b)(^a) 
|ReJA(s)|  ’ 


(6) 


where  Ja(b)  =  dA(B)(—‘^‘^A  +  0).  We  then  replace  the  near-resonant  contributions  Ja(B){s) 
and  Jab{s)  by  constants,  Ja{b)  and  Jab  =  JAB{-i^^A  +  0)  respectively  [see  Eqs.  (4)  and  (5)], 
and  absorb  them  into  the  single-atom  and  two-atom  level  shifts  and  the  respective  decay 
rates.  The  resulting  eigenvalues  and  eigenstates  of  the  singly-excited  system  are  then 


=  ReuJA  ±  ReM, 

lV’5(a)}|{0A})  =  -^{\eA9B)  ±  l5f>ieB))l{0A}), 


(7) 


where  u}a{B)  =  +  Ja(B)-,  M  =  M  +  Jab,  and  s  (a)  denotes  the  symmetric  (antisym¬ 

metric)  eigenstates.  These  two  states  are  split  by  the  RDDI  shift  (which  is  modified  by  the 
structure)  and  give  rise  to  damped  sinusoidal  oscillations  of  the  excitation-transfer  probabil¬ 
ity  from  A  to  5  at  a  rate  of  2ReiV/,  which  varies  as  R~^  at  near-zone  separations,  ojaRIc  <  1 
[2,3]. 

We  shall  be  primarily  concerned  with  the  very  different  behavior  obtained  in  the  case 
of  a  sharp,  narrow  line  or  band  whose  center  of  gravity  is  at  wq.  Consistently,  one  can 
use  in  Eqs.  (4)  and  (5)  the  approximation  Ga{b){<^)  ~  ^ ^'a(B)o^{'^  —  <^o)  and  GabI^^)  ~ 

^  being  the  number  of  near-resonant  modes  in  the  narrow  band  and 
|«.4(S)o|  the  root-mean-square  of  |k.4(b)a|-  This  approximation  holds  if  ^/N\KA(B)o\  is  much 
greater  than  the  band  width.  It  implies  that  we  neglect  all  dissipation^  i.e.,  the  decay 
rates  7^  and  ImM,  originating  from  the  background  density  of  modes  outside  the  line 
(band),  since  they  are  assumed  to  be  much  smaller  than  the  oscillation  frequencies  in  the 
system,  and  we  are  interested  in  times  <  |ImM|"T  As  a  second  simplification,  we 

shall  confine  ourselves  to  the  case  IOa'^aa'^ba  ~  which  holds  for  a  single  near¬ 

resonant  mode  or  a  fully  degenerate  mode,  such  that  ac^a/^ba  is  independent  of  A  (otherwise, 
|Ea«aa«6aI  <  -^|«AO«Bo|)-  We  shall  take  kao  >  0  and  the  sign  of  kbo  to  be  that  of  the 
ratio  kax/kbx-  Following  these  two  simplifications,  the  equation  D{s  =  —iui)  —  0  becomes 


Nk^ 


-  w  - 


AO 


NkI 


ux 


^0, 


U  —  CUB  ~ 


BO 


CU-CUo^ 

+  =  0. 


iO 


—  u;o  y 


(8) 


We  have  now  reduced  the  problem  to  that  of  two  atoms  coupled  to  a  single  mode  via 
coupling  constants  ^/N kaq  and  \/N kbo,  and  to  each  other  via  matrix  element  M.  A  Hamil¬ 
tonian  describing  this  problem  can  be  written  as  a  3  x  3  matrix  in  the  basis  of  the  states 
\^a{B)9b(A),  {0})  and  \9a9b,  lo),  where  0  and  lo  are  the  vacuum  and  single-photon  numbers 
in  the  near-resonant  mode.  We  have  evaluated  in  closed  form  the  eigenvalues  (Fig.  1)  and 
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eigenfunctions  corresponding  to  Eq.  (8),  yet  it  will  be  more  illustrative  to  focus  on  the 
following  distinct  limits: 

(a)  Let  us  assume  that  ~  «so-  This  approximate  equality  of  the  atomic  couplings 
to  the  mode  holds  for  identical  atoms  with  parallel  dipoles  in  the  near  zone  of  separations, 
uaR!  c  In  this  approximation  the  eigenvalues  are  given  by 

Wl,2  «  =  ^(0)5 +a;o  ifi); 

^  —  Wq)^,  (9) 

where  now  Us(^a,)  =  These  results  hold  over  most  of  the  near-zone  separation 

interval  (except  for  the  range  specified  below).  For  |u;s-u;o|  >  '/^{k.ao  +  k.bo),  which  is  the 
case,  e.g.,  for  R  0,  Eq.  (9)  yields  that  and  w_  tend  to  u>s  and  (x^o,  indicating  that  the 

near-resonant  field  mode  is  decoupled  from  the  RDDI-split  symmetric  and  antisymmetric 
states.  Throuhout  the  range  of  validity  of  Eq.  (9),  only  the  antisymmetric-state  eigenvalue 
remains  unperturbed,  since  the  coupling  of  this  state  to  the  mode  2~^''^(kao  —  «:bo)  ~  0 
in  the  near  zone.  By  contrast,  the  symmetric  state  and  the  single-photon  state  become 
hybridized,  giving  rise  to  two  eigenvalues  that  are  split  by  ±f),  the  symmetric-state  vacuum 
Rabi  frequency.  The  above  trends  are  also  reflected  by  the  corresponding  dressed-state 
eigenfunctions 


1-03)  «  |tAa,  {0}),  IV’1,2)  «  iV’±) 

JT  L,.  {0})  1„)] , 

J  i  ±  —  COq  J 


where  c±  —  [1  ±  (cn^  —  u;o)/Q]/2.  We  note  that  I'l/’s)  {0})  and  lffAffB,lo}  coincide  with  |'0+) 
and  for  —  cjqI  >•  V2N(kao  -f  kbo),  but  otherwise  the  symmetric  and  single-photon 
states  are  strongly  mixed  in  |V’±)- 

These  eigenfunctions  can  be  used  to  calculate  the  probability  of  excitation  transfer, 
say,  from  the  initially  excited  atom  A  to  the  initially  unexcited  atom  B,  Pb(^),  and  the 
corresponding  excitation-trapping  probability  PA{t), 


3  2 

PB(A){t)  =  ,  (11) 

t=l 

where  \i^A{B))  —  {0}).  Thus,  the  two-atom  symmetry  determines  the  Rabi 

splitting,  whence  we  find  three  distinct  atomic-state  frequencies,  causing  aperiodic  oscilla¬ 
tions  of  PB{A){t)^  instead  of  sinusoidal  oscillations  in  previously  studied  cases  [2,7],  The 

time  averaged  probabilities,  Pb(a)  =  ZLi\{pB{A)\i’i){'fi\i^A)\  ,  are  approximately  equal, 

•^■4  ^  Rb,  varying  from  3/8  at  —  c<2o|  ■C!  V2N {kao  4-  /cbo)  to  the  free  space  value  1/2  at 
—  wol  »  ^/2N{k,aq  +  «bo)  (note  that  in  this  system  PA{t)  -f-  PB^t)  <  1,  because  of  the 
nonzero  probability  of  the  mode  excitation). 

(b)  Let  us  now  consider  the  effect  of  the  small  near-zone  difference  kao  —  k, bo  which  scales 
approximately  linearly  with  the  separation  R.  This  effect  is  most  remarkable  in  the  range 
of  separations,  corresponding  to  pseudocrossing  (near-equality)  of  two  eigenvalue  curves  in 
Eq.  (9)  (solid  curves  on  Fig.  1),  iox  M  <0  and  023  for  M  >  0.  We  shall 
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assume  below  that  M  is  positive  in  the  near  zone.  The  pseudocrossing  separation  Rc  can  be 
estimated  from  Eq.  (9)  for  the  parameters  that  maximize  ka(B)o-  a  single  mode  confined  to 
the  smallest  possible  volume  ~  (where  A^  =  27rc/w^),  transition  dipole  moment  ~  eo, 
a  being  the  excited  state  radius,  and,  corresppondingly,  ua  ~  e^/(^a).  Using  the  maximal 
estimate  for  RDDI,  M  ~  {eaY  j R?,  this  yields  (i)  R^  ~  yj o.  j {2'k)\ai  where  a  =  1/137, 

if  \ujo-<jja\  ^  V''^(k.4o+«bo);  (ii)  Rc  ~  \ioJo-uJA  >  \/^(«xo+«Bo)- 

In  either  case,  Rc  A^,  i.e.,  the  pseudocrossing  occurs  well  within  the  near  zone.  Note  that 
for  certain  dipole  orientations  M  is  further  reduced  and  Rc  can  become  even  smaller. 

Drastic  modifications  of  the  interatomic  coupling  occur  near  Rc  (as  compared  to  the 
cases  discussed  above),  due  to  the  strong  competition  of  RDDI  and  Rabi  splittings.  The 
eigenvalues  in  Eq.  (9)  are  then  replaced  by  the  more  accurate  solutions  of  Eq.  (8) 


ui  «  a;+,  u;2,3  ~  -(^<^-  +  Wa  ±  Cl'), 


(12) 


Here  |Vo|  equals  the  minimal  splitting  io-i  —  and  also  determines  the  width  of  the  pseu¬ 
docrossing  interval,  which  is  defined  as  follows:  ARc  =  |i?i  —  R2I5  where  ijJa-{Ri,2)  =  |Eoi 
and  Ufa-  =  Iwa  —  w_|.  For  two  atoms  far  from  a  node  of  a  sinusoidal  mode,  Vo  ~ 
N kao\i^ao  —  i^bo\/{V8Nkao  +  |<^o  “^^sl)  3-nd  is  proportional  to  Rc  at  near-zone  separations 
[11]. 

The  eigenfunction  \tp+)  is  not  affected  by  the  pseudocrossing,  whereas  |^_)  and  |^a)  are 
strongly  mixed  near  Rc,  forming 


_ ^ 

±  (uj- 


#a,{0}) 


(13) 


where  62,3  =  [l±(u;_— a;a)/fi']/2.  This  mixing  implies  the  complete  breaking  of  the  symmetry 
[Eq.  (7)],  which  characterizes  the  two-atom  system  subject  to  RDDI  in  the  weak-coupling 
case.  At  R  —  Rc,  IV’2,3)  =  2“^/^(|V’_)±|V’a,  {0})).  For  sufficiently  large  and  positive  detuning 
Wo  -Ws  ~  D,  c+  ^  0,  yielding  [see  Eq.  (10)]  [^’2)  \^a9b,{^]),  and  \i)^z)  {0}), 

which  means  that  the  excited  states  of  atoms  A  and  B  become  uncoupled,  due  to  interference 
of  \ips)  and  |0a)-  The  corresponding  time-averaged  value  of  the  excitation  transfer  probability 
exhibits  strong  suppression  in  the  pseudocrossing  interval,  Pb{R  —  Rc)  =  3c^/8  1. 

Concurrently,  Pa{R  =  Rc)  =  1  —  c+  +  (3/8)c^,  tending  to  1  with  the  increase  of  wq  —  Wj 
(solid  curves  in  Fig.  2).  This  means  that  the  excitation  is  then  strongly  trapped  at  the  initial 
atom,  which  reflects  the  decoupling  of  the  atomic  excited  states.  These  results  sharply 
contrast  those  outside  the  crossing  region  [limit  (a)  above]. 

(c)  Finally,  let  us  examine  the  limit  of  large  differences  between  the  couplings:  k*  = 
—  «bo|  ~  +  «^Bo|-  For  atoms  with  non-parallel  dipoles,  this  limit 

can  occur  at  any  separation.  The  analysis  of  Eq.  (8)  in  this  limit  shows  that  the  results  of 
limit  (b)  apply  to  the  pseudocrossing  of  Wa  with  a;_,  if  uq  —  u a  >■  «:s(a)-  The  only  salient 
difference  with  limit  (b)  is  in  the  width  of  the  pseudocrossing  region.  According  to  the 
discussion  following  Eq.  (12)  this  region  is  defined  by  Wa_(i?)  =  |wo(R)  —  a;_(R)|  ^  ]Vo(i?)]. 
Under  the  condition  wq  —  ua  ~  Ka,  we  find  that  uJa-{R)  increases  with  R  for  R  >  Rc, 


until  it  attains  its  upper  bound  Nk1I{ujq  —  o;^).  On  the  other  hand,  |Vo|  decreases  for 
R  >  Rc  down  to  its  lower  bound  NksK^Kojo  —  (actually  this  bound  is  slowly  varying 
due  to  the  spatial  dependence  of  but  this  does  not  change  |ho|  appreciably).  Since  these 
two  bounds  coincide,  the  pseudocrossing  region  can  extend  from  i?  <  to  i?  > 

(dashed  curves  in  Fig.  1)  (!).  Throughout  this  region,  Eqs.  (12)  and  (13)  hold  fully,  i.e., 
excitation  transfer  is  strongly  suppressed  by  interference  of  and  |^a)  (dashed  curves  in 
Fig.  2).  Similar  results  hold  also  for  cu.4  -  wq  >  Ks(a)- 

The  analytical  treatment  outlined  above  is  in  full  agreement  with  the  numerical  analysis 
of  the  exact  Eq.  (3)  (as  shown  in  Eigs.  1  and  2).  Both  treatments  demonstrate  that  two 
identical  atoms  interacting  via  a  near-resonant  narrow-linewidth  mode  (or  degenerate  band) 
and  an  off-resonant  reservoir  can  exhibit  a  much  richer  variety  of  spectral  and  dynamical 
features  than  what  is  currently  known.  It  must  be  viewed  as  a  system  of  three  mutually 
coupled  excited  states  (as  opposed  to  two  such  states  in  previously  studied  models  [2,7]). 
The  dipole-dipole  (RDDI)  coupling  and  the  mode-induced  Rabi  splitting  are  inseparable  in 
this  system.  The  interatom  coupling  results  from  mixing  of  all  three  states,  which  is  most 
striking  for  separations  where  the  two  predominantly-populated  levels  nearly  cross.  Then 
the  competing  RDDI  and  Rabi  splittings  cause  strong  interference  of  the  symmetric  and 
antisymmetric  two-atom  excited  states,  leading  to  decoupling  of  single-atom  excited  states. 
This  occurs  at  near-zone  (quasimolecular)  separations  (as  small  as  lO'^A^)  and  corresponds 
to  suppression  of  interatom  excitation  transfer,  compared  to  the  electrostatic  near-zone  limit 
of  RDDI.  ■ 

The  present  predictions,  particularly  the  suppression  of  excitation  transfer  at  quasimolec¬ 
ular  separations,  may  be  important  in  various  systems  of  two  atoms,  molecules,  or  excitons 
within  high-(5  resonators.  Dielectric  sphere  surfaces  and  Bragg  resonators  are  particularly 
suitable  configurations.  Control  of  the  initial  state  of  each  partner,  even  at  short  separa¬ 
tions,  can  be  achieved  by:  (i)  excitation  of  only  one  partner,  say  A,  to  [e^),  via  another 
state  |n,4)  which  differs  horn  its  counterpart  in  B,  although  the  transitions  je^)  — >  |^^)  and 
1^5)  ^  are  identical  (this  option  is  common  in  molecules  or  impurity- bound  excitons); 

(ii)  single-partner  excitation  localized  on  a  10^  nm  scale  (by  electronic,  XUV  or  near-held 
optical  techniques);  (iii)  slow  (ultracold)  atom  collisions  in  a  cavity,  with  one  atom  initially 
excited  outside  the  cavity.  The  dipole  orientations  are  controllable  by  the  initial  excitation 
and  (or)  by  magnetic  splitting  of  Zeeman  sublevels. 
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coupled  to  a  near-resonant  mode.  Solid  curves  -  exact  eigenvalues  as  a  function  of  the  normalized 
separation  RJXa  for  N  =  1,  M  =  (7^/2)(cx;^i?/c)~^,  where  7^  is  the  free  space  decay  rate  of  the 
atom  excited  state,  uq  —  oja  =  10^7^,  ka(B)  =  10^7a  sin(w^r^(B)/c),  and  =  O.SA^.  Note  the 
narrow  pseudocrossing  region.  Dashed  curves  -  idem,  for  ka  =  IOOO7A  and  kb  =  5007^4.  The 
pseudocrossing  region  is  now  very  broad. 
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as  in  Fig.  1:  solid  curves  -  for  narrow  interval  of  excitation-transfer  suppression;  dashed  curves  - 
for  broad  interval  of  suppression. 
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Appendix  D: 

Self-Induced  Transparency  in  Bragg  Reflectors: 

Gap  Solitons  Near  Absorption  Resonances 

Alexander  Kozhekin  and  Gershon  Kurizki 
Chemical  Physics  Department,  Weizmann  Institute  of  Science,  Rehovot  76100,  Israel 

We  show  that  pulse  transmission  through  near-resonant  media  embedded 
within  periodic  dielectric  structures  can  produce  self-induced  transparency 
(SIT)  in  the  band  gap  of  such  structures.  This  SIT  constitutes  a  principally 
new  type  of  gap  soliton. 

42.50Rh,  42.25.Bs,  03.40.Kf,  78.66.-w. 

Self-induced  transparency  (SIT),  namely,  solitary  propagation  of  electromagnetic  (EM) 
pulse  in  near-resonant  media,  irrespective  of  the  carrier-frequency  detuning  from  resonance, 
is  one  of  the  most  striking  and  important  effects  of  nonlinear  optics  [1,2].  It  reflects  the 
essence  of  driven  two-level  atom  dynamics,  which  is  described,  in  the  soliton  frame,  by  a 
pendulum  equation  for  the  pulse  area  6  (the  sine-Gordon  equation).  If  the  pulse  duration 
is  much  shorter  than  the  transition  (spontaneous-decay)  lifetime  (Ti)  and  dephasing  time 
(r2),  and  0  is  a  multiple  of  27r,  then  pulse-area  conservation  gives  rise  to  SIT,  corresponding 
to  fe-emission  of  the  absorbed  radiation  in-phase  with  the  driving  field. 

One  of  the  standard  tacit  requirements  for  SIT  is  uniformity  of  the  medium.  Indeed, 
one  would  expect  that  partial  reflection  of  the  field  in  a  non-uniform,  e.g.,  layered,  medium 
should  destroy  SIT,  because  the  pulse  area  is  then  split  between  the  forward  and  backward 
(reflected)  waves,  and  is  no  longer  conserved  for  each  wave.  This  expectation  seems  to  be 
supported  by  treatments  of  a  single  thin  resonant  film  [3]  or  periodic  array  of  such  films  [4], 
which  yield  bistable  (two- valued)  transmission  of  the  incident  pulse,  owing  to  the  coupling  of 
the  forward  and  backward  waves.  Should  we  then  anticipate  severely  hampered  transmission 
through  a  medium  whose  resonance  lies  in  a  reflective  spectral  domain  (photonic  band  gap) 
of  a  periodically-layered  structure  (a  Bragg  reflector)?  In  this  Letter  we  show  that,  contrary 
to  such  expectations,  it  is  possible  for  the  pulse  to  overcome  the  band-gap  reflection  and 
produce  SIT  in  a  near-resonant  medium  embedded  in  a  Bragg  reflector. 

The  predicted  SIT  propagation  is  a  principally  new  type  of  gap  soliton,  which  does 
not  obey  any  of  the  familiar  soliton  equations,  such  as  the  non-linear  Schrodinger  equation 
(NLSE)  or  the  sine-Gordon  equation.  Its  spatio-temporal  form,  intensity  dependence  and 
transmission  mechanism  are  shown  here  to  be  quite  unique.  Nevertheless,  it  shares  several 
common  features  with  the  extensively  -  studied  gap  solitons  [5-8]  or  with  ultrashort  pulses 
[9]  in  Kerr-nonlinear  Bragg  reflectors. 

We  consider  the  propagation  of  an  EM  pulse  through  a  medium  consisting  of  two-level 
systems  (TLS)  embedded  within  a  one-dimensional  periodic  dielectric  structure,  e.g.,  a  multi¬ 
layered  dielectric  mirror.  Our  starting  point  is  the  Maxwell  equation  for  the  field  E 

dfi  '  ' 

driven  by  the  current  of  the  total  polarization  Ptot  “  ^in  +  The  linear  part  of  the 
polarization  Pj-jj  =  Xlin^  characterized  by  a  linear  refractive  index  =  (1  + 
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fundamental  period  d,  =  cq  +  Em=i  cos{2mkz),  where  k  =  %  j d  and 
Acm  is  the  variation  of  m-th  harmonic  of  the  dielectric  index.  The  non-linear  polarization 
is  the  near-resonant  response  of  the  TLS. 

We  may,  analogously  to  the  theory  of  distributed  feedback  lasers  [10-12],  decompose  the 
total  field  and  non-linear  polarization  into  forward  (F)  and  back-reflected  (B)  components 


E{z,t)  = 
^nl(^’0  = 


VF{z,t)e^'^^  +  VB{z,t)e-^'‘^ 


^-iWgct 


(2a) 

(2b) 


assuming  that  the  carrier  frequency  u  lies  near  the  center  of  the  lowest  fundamental  band 
gap  Ugc  =  kclno^  where  no  =  The  small  detuning  from  the  gap  center  \ugc  —  u;|  <C  cj 
will  be  considered  as  phase  modulation  of  the  complex  amplitude  £f(B)-  Under  the  weak- 
reflection  assumption  (|Aem|  Cq)  we  may  drop  spatially  fast- varying  components  (varying 
on  the  scale  of  a  wavelength)  of  E  and  ^nl>  and,  consistently,  m  >  1  terms  of  We 

then  obtain  the  following  coupled-mode  equations  for  the  Rabi  frequencies  corresponding  to 
the  slow-varying  field  amplitudes,  £If[B)  =  fJ-£F{B)lk,  where  n  is  dipole  moment  of  the  TLS 
transition  [12] 


\  no  dz  dt ) 


f^F(B)  —  icK/noClB{F)  +  R  '^E’f{B) 


(3) 


Here  and  hereafter  the  upper  (lower)  sign  corresponds  to  the  first  (second)  subscript.  The 
first  term  on  the  r.h.s.  of  (3)  describes  the  forward-to-backward  wave  coupling  via  Bragg 
reflection  with  characteristic  reflection  (attenuation)  length  l//c,  where  k  =  ^Aei/4eo.  The 
second  term  is  proportional  to  the  TLS  polarization  current,  and  cr) / {nlU) 

is  the  inverse  square  of  the  cooperative  resonant  absorption  time  r^,  cr  being  the  density  of 
the  TLS.  Note  that  k  is  positive  or  negative  for  TLS  embedded  in  regions  with  the  higher 
or  lower  refractive  index  (Aei  >  0  or  Aci  <  0),  respectively. 

In  treatments  of  bidirectional  field  propagation  in  media  with  arbitrary  spatial  distribu¬ 
tion  of  near-resonant  atoms  [13,14],  the  Bloch  equations  for  the  population  inversion  w{z,t) 
and  polarization  are  entangled  in  a  fashion  which  leads  to  an  infinite  hierarchy  of  equa¬ 
tions  for  successive  spatial  harmonics.  The  truncation  of  this  hierarchy  can  only  be  justified 
by  phenomenological  arguments.  Here  we  avoid  this  complication  by  restricting  the  near¬ 
resonant  species  distribution  to  thin  layers  (much  thinner  than  the  resonant  wavelength) 
with  the  same  periodicity  d  as  the  dielectric  structure.  Then  the  nonlinear  polarization 
envelopes  may  be  decomposed  as  [4] 

=  )£’=■*'  (4) 

j 


Here  P^^{zj,t)  is  the  nonlinear  polarization  of  the  TLS  in  the  j-th  layer  [zj  =  jd).  The 
Bloch  equations  then  assume  the  form 

dtP^^{zj,  t)  =  w{Zg,  t) 

-  i(cjgc  -  ui2)P^i(zj,  t)  (5a) 

dtw(zj,t)  =  -lp*^(zg,i) 

+  c.c.  (5b) 
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where  to  is  the  population,  inversion  and.  u)gc  —  ^12  is  the  detuning  of  the  TLS  resonance 
frequency  U12  from  the  gap  center.  The  periodicity  of  the  TLS  positions  which  satisfy  the 
Bragg  condition,  i.e.,  exp  {±2ikzj)  =  1  yields  Vf  =  Vb='P  in  (4).  Upon  summing  over  all 
layers  j,  the  Bloch  equations  (5)  can  be  written  under  the  Bragg  condition  in  the  following 
closed  form,  without  resorting  to  harmonic  expansion 

dtV{z^  t)  =  w{z,  t)  {(If  +  (Ib)  —  —  '^\2)'P{z,  t)  (6a) 

dtw{z,t)  = --\P*{z,t){(lF +  ^3)  c.c.  ]  (6b) 


In  an  attempt  to  further  simplify  the  Maxwell-Bloch  equations,  we  first  convert  the 
variables  to  the  dimensionless  form  t  =  t/xc  ;  (  =  nozj{cTc)  ;  rj  =  kTc  and  8  =  {ugc  —  n>i2)'>'c- 
We  now  rewrite  Eqs.(3)  and  (6)  for  the  sum  and  difference  of  the  forward  and  backward  field 
envelopes  =  Tc{0.f  +  (^b)  and  S_  =  r<,(U/r  -  Ub),  and  obtain  by  simple  manipulations 


S+  =  2drV  +  t2rjV  - 
S_  =  -2dcV  -  77^S_ 


(7a) 

(7b) 


Although  and  S_  now  obey  separate  equations,  they  are  still  coupled  via  V ,  which 
satisfies  the  Bloch  equations 


drV  =  -  iSV  (8a) 

d,w  =  -^{v*^++VJ:i)  (8b) 

We  emphasize  again  the  crucial  role  of  the  assumption  that  the  TLS  layers  are  much  thinner 
than  a  wavelength  and  satisfy  the  Bragg  condition.  Without  this  assumption  we  could  not 
have  obtained  (7)  and  (8),  which  are  closed  in  V  and  w  (in  contrast  to  Ref.  [13,14]). 

Equations  (7)  and  (8)  cannot  be  reduced  to  any  familiar  soliton  equation.  Our  main  idea 
is  to  try  for  the  above  equations  a  phase-modulated  27r-soliton  SIT  solution 


2+ 


exp  (i(a(C  —  Ar)) 
°cosh  [(^(C/'w  —  t)] 


(9) 


where  A  =  {u  —  ujgc)Tc^  Ao  is  the  amplitude  of  the  solitary  pulse,  ^  its  width  and  u  its  group 
velocity  (normalized  to  c). 

Substituting  drV  from  Eq.(8a)  into  Eq.(7a),  we  may  express  V  in  terms  of  E+  and 
the  population  inversion  w.  Then,  upon  eliminating  V  and  using  Eq.(9),  we  can  integrate 
Eq.(8b)  for  the  population  inversion  w,  obtaining 


Ao^(A  —  a/u) _ 1 _ 

2{8  —  k)  cosh^  [/3(C/«  -  t)] 


(10) 


Using  these  explicit  expressions  for  V  and  w  in  (7a)  and  (8a),  we  reduce  our  system  to  a  set 
of  algebraic  equations  for  the  coefficients  a,  A  that  determine  the  spatial  and  temporal  phase 
modulation,  and  the  pulse  width  jd  as  functions  of  the  velocity  u.  The  soliton  amplitude  is 
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then  found  to  satisfy  |y4o|  =  2/9,  exactly  as  in  the  case  of  usual  SIT  [2].-  This  implies,  by 
means  of  Eq.(9),  that  the  area  under  the  E_|-  envelope  is  27r. 

Let  us  consider  the  most  illustrative  case,  when  the  atomic  resonance  is  exactly  at  the 
center  of  the  optical  gap,  5  =  0.  Then  the  solutions  for  the  above  parameters  are  [15] 


a 


7/  1  —  Zv? 

2u  1  — 


A  = 


7?1 


2  1 


u 


=  \Ao\Vi  = 


_  8u^(l  —  u'^)  —  7/^(1  + 
4(1 


u 


2\2 


(Ua) 

(11b) 


In  the  frame  moving  with  the  group  velocity  of  the  pulse,  C,'  =  (  -  ut,  the  temporal 
phase  modulation  will  be  [au  —  A)r,  which  is  found  from  Eq.(ll)  to  be  equal  to  —rjT.  Since 
7]  =  KTc  is  the  (dimensionless)  gap  width,  this  means  that  the  frequency  is  detuned  in  the 
moving  frame  exactly  to  the  band-gap  edge.  The  band-gap  edge  corresponds  (by  definition) 
to  a  standing  wave,  whence  this  result  demonstrates  that  such  a  pulse  is  indeed  a  soliton, 
which  does  not  disperse  in  its  group-velocity  frame. 

The  allowed  range  of  the  solitary  group  velocities  may  be  determined  from  Eq.(llb) 
through  the  condition  >  0  for  a  given  t/,  as  illustrated  in  Fig.l.  The  same  condition 
implies  |7;|  <  77  max,  where 


2  _  8w^(l  —  u^) 

’  ^  (1  +tl2p 


(12) 


It  follows  from  (12)  that  the  condition  for  SIT  is  [77 1  <  1,  77  max  =  1  corresponding  to 
u  =  l/\/3.  This  condition  means  that  the  cooperative  absorption  length  crdno  should  be 
shorter  than  the  reflection  (attenuation)  length  in  the  gap  1/ac,  i.e.,  that  the  incident  light 
should  be  absorbed  by  the  TLS  before  it  is  reflected  by  the  Bragg  structure.  In  addition, 
both  these  lengths  should  be  much  longer  than  the  light  wavelength  for  the  weak- reflect  ion 
and  slow-varying  approximation  to  be  valid. 

From  Eq.(7b)  we  find 


S_ 


(13) 


and  the  equation  for  population  inversion,  obtained  from  Eq.(lO) 


tu  =  —  1  -f- 


cosh^  [/3{C/u  —  r)] 


(14) 


The  envelopes  of  both  waves  (forward  and  backward)  propagate  in  the  same  direction; 
therefore  the  group-velocity  of  the  backward  wave  is  in  the  direction  opposite  to  its  phase- 
velocity!  This  is  analogous  to  climbing  a  descending  escalator. 

Analogously  to  Kerr- nonlinear  gap  solitons  [5,6],  the  real  part  of  the  nonlinear  polar¬ 
ization  RePjj  creates  a  traveling  ’’defect”  in  the  periodic  Bragg  reflector  structure  which 
allows  the  propagation  at  band-gap  frequencies.  The  real  part  of  the  nonlinear  polariza¬ 
tion  is  governed  by  the  frequency  detuning  from  the  TLS  resonance.  Exactly  on  resonance 
(which  we  here  take  to  coincide  with  the  gap  center)  A  =  5  =  0,  =  0,  and  qur 
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solutions  (11)  yield  imaginary  values  of  the  velocity  u  and  modulation  coefficient  a.  The 
forward  field  envelope  then  decays  with  the  same  k  exponent  as  in  the  absence  of  TLS  in  the 
structure.  Because  of  this  mechanism,  SIT  exists  only  on  one  side  of  the  band-gap  center, 
depending  on  the  sign  of  k  in  Eqs.(3),  i.e.,  on  whether  the  TLS  are  in  the  region  of  the 
higher  or  the  lower  linear  refractive  index.  This  result  may  be  understood  as  the  addition  of 
a  near-resonant  non-linear  ”refractive  index”  to  the  modulated  index  of  refraction  of  the  gap 
structure.  When  this  addition  compensates  the  linear  modulation,  then  soliton  propagation 
is  possible  (Fig. 2).  On  the  "wrong”  side  of  the  band-gap  center,  soliton  propagation  is  for¬ 
bidden  even  in  the  allowed  zone  because  the  nonlinear  polarization  then  cannot  compensate 
even  for  a  very  weak  loss  of  the  forward  field  due  to  reflection. 

The  soliton  amplitude  and  velocity  dependence  on  frequency  detuning  from  the  gap 
center  (which  coincides  with  atomic  resonance)  are  illustrated  in  Fig.l.  They  demon¬ 
strate  that  forward  soliton  propagation  is  allowed  well  within  the  gap,  for  A  satisfying 
(1  —  -y/l  —  r]‘^)lrj  <  A  <  (1  -|-  a/1  —  In  addition  to  frequency  detuning  from  reso¬ 

nance,  the  near-resonant  gap  soliton  possesses  another  unique  feature:  spatial  self-phase 
modulation  aC  of  both  the  forward  and  backward  field  components. 

To  check  our  analytical  solutions  (9)-(13)  we  have  compared  them  with  numerical  simula¬ 
tions  of  Eq.(3)  using  the  numerical  method  developed  in  [4].  As  the  launching  condition,  we 
take  the  incident  wave  in  the  form  Ep  =  A  exp  [i(a;  —  <jOgc){t  —  to)]/  cosh  [/3(t  —  to)/'^c]  with¬ 
out,  a  backward  wave  {£b  =  0)  at  ^  =  0.  By  varying  the  detuning  to  —  Ugc  and  amplitude  A 
we  investigate  the  field  evolution  inside  the  structure.  When  these  parameters  are  close  to 
those  allowed  by  Eqs.(ll),(12),  we  observe  the  formation  and  lossless  propagation  of  both 
forward  and  backward  soli  ton-like  pulses  with  amplitude  ratios  predicted  by  our  solutions 
(Fig. 3. a).  By  contrast,  exponential  decay  of  the  forward  pulse  in  the  gap  is  numerically 
obtained  in  the  absence  of  TLS  (Fig.3.b). 

The  observation  of  the  predicted  SIT  at  band-gap  frequencies  requires  high  values  of 
the  dipole  moment  p  and  high  density  of  the  TLS,  in  order  to  achieve  short  Tc  along  with 
large  T2  and  T\  times  (to  avoid  dephasing  and  energy  losses  by  incoherent  processes).  The 
most  adequate  system  for  experimental  observation  of  this  effect  appear  to  be  excitons  in 
semiconductors.  Let  us  consider  a  periodic  array  of  12-nm-thick  GaAs  quantum  wells  (A  = 
806nm)  separated  by  A/2  non-resonant  AlGaAs  layers  [16].  In  this  system  p  ~  10~^®cgse 
and  area  density  concentration  cr  ~  10*  cm“^  are  achievable  (corresponding  to  an  average 
bulk  density  of  ~  10^*  cm~*)  which  yields  ~  10“^*s.  The  relaxation  time  (at  2°K)  is 
then  T2  ~  lOps  [16].  A  solitary  pulse  duration  of  the  order  of  ^  Ips,  i.e.,  much  shorter  then 
the  dephasing  time  T2  corresponds  to  a  pulse  width  ,5  ^  0.1.  From  Eq.(ll)  we  find  that 
SIT  in  this  structure  requires  that  7/  <  0.99,  corresponding  to  a  band-gap  reflection  length 
1/k  >  100  A,  i.e.,  Aci/eo  ^  0.01.  An  alternative  choice  may  be  bound  I2  excitons  in  CdS, 
which  yield  similar  Tc  and  T2  at  lower  exciton  densities,  (controlled  by  the  donor  impurity 
concentration)  and  are  therefore  more  appropriable  for  the  TLS  description. 

To  conclude,  we  have  demonstrated  the  possibility  of  solitary  pulse  transmission  through 
a  one-dimensional  band-gap  region,  by  means  of  near-resonant  polarization  effects.  In  com¬ 
parison  with  Kerr-nonlinear  gap  solitons,  near-resonant  gap  solitons  have  several  unique 
features:  temporal  and  spatial  pha^e  modulation,  the  detuning  conditions  and  the  result¬ 
ing  velocity  and  amplitude  threshold.  Their  salient  advantage  is  stability  with  respect  to 
absorption.  By  contrast,  strong  absorption  is  a  severe  problem  associated  with  a  large  non- 
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linear  Kerr  coefficient.  As  regards  applications,  the  sensitivity  to  launching  conditions  on 
the  carrier  frequency  and  pulse  shape  can  make  the  predicted  gap  SIT  an  effective  filter  for 
signal  transmission. 
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and  the  Minerva  Foundation,  Germany. 
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FIG.  1.  Dependence  of  the  solitary  pulse  velocity  (solid  line)  and  amplitude  (dashed  line)  on 
frequency  detuning  from  the  gap  center  for  tj  =  0.7.  At  the  gap  edge  (dotted  line)  u  =  l/v/S  and 


curve)  in  a  structure  of  periodically  alternating  layers.  This  modulation  can  be  canceled  by  the 
near-resonant  nonlinear  response  Ree^^  (inset),  if  it  has  the  opposite  sign  to  Aci  at  the  TLS 
positions. 


FIG.  3.  Numerical  simulations  of  the  intensities  of  ’’forward”  waves  in  the  gap:  (a)  when 

Eqs.(ll),  (12)  are  obeyed  (r?  =  0.7,  group  velocity  u  ~  0.3);  (b)  without  TLS  (same  77  and  incident 
pulse  as  in  (a)). 
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Appendix  E: 

Optical  ’’Multi-Excitons”:  Quantum  Gap  Solitons  in  Nonlinear 

Bragg  Reflectors 

Ze  Cheng  and  Gershon  Kurizki 

Department  of  Chemical  Physics,  Weizmann  Institute  of  Science,  Rehovot  76100,  Israel 

We  have  accomplished  the  Bethe-ansatz  solution  for  pairwise  interacting 
quanta  within  the  effective-mass  regime  of  band-gap  propagation  in  nonlinear 
Bragg  reflectors.  Our  theory  predicts  a  new  kind  of  collective  excitations  of 
the  electromagnetic  field  dressed  by  such  media,  namely,  optical  multi-exciton 
(OME)  complexes  (or  condensates),  which  are  the  quantum  states  associated 
with  gap  solitary  waves.  Their  existence  should  be  manifested  by  the  discrete 
spectrum  of  band-gap  transmission  as  a  function  of  the  incident  photon  num¬ 
ber  and  by  the  multi-exponential  faU-off  of  intensity-intensity  correlations  on 
a  0.1  mm  scale.  OMEs  should  have  advantageous  stability  properties. 


Quantum  effects  of  light  propagation  in  nonlinear  fibers  have  been  the  subject  of  exten¬ 
sive  studies  in  recent  years  [1].  These  studies  have  established  that  optical  solitons  which 
are  classically  described  by  the  nonlinear  Schrodinger  equation  (NLSE)  [2],  have  quantum 
analogs  in  the  form  of  superposition  of  mutually-bound  (spatially-correlated)  multiphoton 
states.  The  description  of  such  bound  multiphoton  states  is  given  by  the  Bethe  ansatz 
solution  of  the  second-quantized  NLSE  [3]. 

Here  we  address  the  hitherto  unexplored  quantum  regime  of  optical  propagation  in  for¬ 
bidden  spectral  bands  (band-gaps)  of  Kerr-nonlinear  one-dimensional  (ID)  Bragg  reflectors, 
wherein  transmission  of  solitary  waves  (gap  solitons)  [4]  and  ultrashort  pulses  [5]  has  been 
predicted  classically.  By  extending  the  Bethe  ansatz  solution  to  such  systems,  we  reveal 
and  investigate  the  striking  analogy  between  mutually-bound  multiphoton  states  of  quan¬ 
tum  gap  solitons  and  multi-exciton  complexes  (excitonic  molecules)  in  semiconductors  [6]. 
These  optical  multi-excitons  (OME)  can  be  formed  by  the  combined  effect  of  (i)  photon 
effective  masses,  which  are  endowed  by  the  periodic-structure  dispersion;  (ii)  Kerr-nonlinear 
interband  photon  attraction,  and  (iii)  their  intraband  attraction  or  repulsion,  also  caused  by 
Kerr  nonlinearity.  The  main  manifestations  of  OMEs  are  predicted  here  to  be:  (1)  a  series  of 
discrete  transmission  lines  in  the  optical  band  gap  of  the  Bragg  reflector,  which  correspond 
to  increasing  numbers  of  photon  pairs  and  range  from  a  single  exciton  to  the  ’’ionization” 
threshold;  (2)  multi-exponential  fall-off  of  two-point  intensity-intensity  correlations.  Both 
features  exhibit  dramatic,  unparalleled  dependence  on  the  refractive-index  modulation  in 
the  Bragg  reflector. 

The  effective  Hamiltonian  in  terms  of  the  magnetic  and  electric  field  operators  B  and  E 
in  a  ID  periodic  Kerr-nonlinear  dielectric  structure  can  be  written  as 


+  ^-^n\z)E^E  -f 

2/xo  2  ^  ^  4xp2 


(1) 


where  :  :  denotes  normal  ordering.  The  Hamiltonian  in  Eq.(l)  is  similar  to  that  used  in 

studies  of  quantum  solitons  in  nonlinear  fibers  [7]  except  that  the  dielectric  index  7z^(z)  and 
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the  nonlinear  susceptibility  are  now  periodic  in  Here  tt/?^  is  the  interaction  region 

area,  determined  by  the  optical  beam  waist.  On  making  the  standard  assumption  that  only 
the  two  bands  bordering  a  certain  band  gap  are  involved  in  the  nonlinear  process  [6],  we 
can  quantize  the  Hamiltonian  in  Eq.(l)  in  the  two-band  Bloch-function  basis,  denoting  the 
lower  ('Valence”)  and  upper  (’’conduction”)  bands  by  -f  and  respectively.  The  quantized 
Hamiltonian  reads 


a—± 


^^cxk^al^Q'ak  +  L  ^  Vaq{k^  k 

k  k,k’ ,q 


k,k'  ,q 


k'-\-q  ^-\-,~k^^~{-,~k^-,k-{-q  • 


(2) 


Here  L  is  the  structure  length,  Vaq{k^k')  are  the  intraband  matrix  elements  of  the  Kerr- 
nonlinear  ’’potential”,  and  Viq{k,k')  are  the  corresponding  interband  matrix  elements. 

We  are  interested  in  photon  properties  near  the  extrema  of  the  two  bands,  where  k  =  ko. 
We  can  therefore  use  the  well-known  effective-mass  approximation  [6] 


—  hu)-^ 


n%k  -  ko)^ 
2m± 


(3) 


where  m±  =  T^/{d‘^oj±k/ are  the  photon  effective  masses  in  the  two  bands.  In  a 
structure  of  alternating  layers  with  refractive  indices  ni,  no  and  thicknesses  ax,  a2,  in  the 
case  nxai  =  n2CL2  and  for  the  second  lowest  band  gap,  the  photon  effective  masses  are  found 
to  satisfy 


^TT^i  CLi 

c(ax  T  0,2^"^ 


fl  ,  1 

{n2 

,  ni\ 

1  +  “ 
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+  — 

L  2 

Vnx 

n2y  J 

sm 


2niaiu;± 


(4) 


One  sees  that  |m_  —  m^\  scales  with  the  band  gap  width  when  \n2  —  ni  |  <C  ni^2-  Consistently 
with  the  effective-mass  approximation,  we  replace  the  potential  matrix  elements  in  Eq.(2) 
by  their  values  V±  and  W  at  ^  =  0  and  k  =  k'  =  ko.  As  we  shall  be  working  in  the  coordinate 
representation,  we  take  the  Fourier  transform  of  the  annihilation  operators  dak,  obtaining 
the  following  position- dependent  field  operators 


<^±(z)  =  ii-C2Ea±,e‘("-^°)F  (5) 

k 

The  Hamiltonian  (2)  can  now  be  rewritten  as 

^  =  ^(a;+-fa;_)iV  +  ^„  (6a) 

where  A  =  /  (f>l^(f>adz  is  the  single-band  photon-number  operator,  and 

H,  =  JdzlYl  -b  V^4>UUa<i>a]  +  V,ct>Ul<t>^4>\  .  (6b) 

V  a=±  J 

The  ’’kinetic  energy”  terms  depend  on  the  structure  dispersion,  whereas  the  ’’potential 
energy”  signs  are  those  of  the  Kerr  nonlinearity  i.e.,  positive  (negative)  for  self- focusing 
(self-defocusing)  media. 


42 


Since  Hs  commutes  with  N,  we  seek  simultaneous  eigenstates  of  both  operators.  Such 
eigenstates  can  be  obtained  exactly  by  the  Bethe  ansatz  method,  by  extending  Thacker’s 
integrable  quantum- field  theory  [8]  to  the  present  two-band  system.  Since  each  band  contains 
an  equal  arbitrary  photon  number  N,  the  Hamiltonian  in  Eqs.(6)  admits  only  even  numbers 
of  quanta.  We  first  construct  the  exact  four-quantum  eigenstate,  and  then  generalize  it  to 
arbitrary  even  quantum  numbers  2,N.  The  resulting  2A^-quantum  eigenstate  is  obtained 
from  a  tedious  calculation  according  to  the  aforementioned  procedure.  It  is  given  in  its 
unnormalized  form,  by 


1$2V  >—  /  •••? -^/V±)  n  JX  \4>ai^ja)d'^ja  |0  >5 

j=l  a=± 


(7) 


where  the  spatial  envelope  (wave  function)  is 
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(8) 


The  wave  function  'I'(2i±, 2a/±)  is  symmetric  with  respect  to  the  exchange  of  any  two 
subscripts  ]  and  /  for  photonic  coordinates,  consistently  with  the  bosonic  nature  of  photons. 
The  application  of  the  Hamiltonian  H  to  the  eigenstate  |$2Af  >  yields  the  eigenenergy 


N  (  \ 

E2V  =  Nh{LJ+  +  u;_)  +  ^  f  r,  j  • 


(9) 


In  the  above,  all  the  wave  vectors  are  measured  relative  to  the  band-edge  wave  vector  ko. 

We  are  looking  for  bound-state  solutions  of  (8),  wherein  the  exponential  factors  are  real, 
and  fall  off  with  the  separation  of  photons  from  different  bands  \zj^  —  Zj_|  or  from  the  same 
band,  \zja  -  It  can  be  checked  that  the  wave  vectors  kj^  in  (8)  must  then  satisfy 


kj±  — 


m± 


-K±i 


—  m-f 


(m_  —  m^)h 


,  =F  i^(iV  -  2i  +  1), 


(10) 


where  a  =  ±  is  chosen  according  to  whether  the  Kerr  nonlinearity  is  positive  (self-focusing) 
or  negative  (self-defocusing),  and  NK  is  the  total  composite  momentum,  which  can  be 
expressed  as 


N 


NK  =  Y.{l=,+  +  <=!-)■ 

j=l 


(11) 


Equation  (11)  indicates  that  bound  states  consist  of  N  ’’excitons”,  i.e.,  pairs  of  ’’conduction”- 
and  "valence” -band  photons.  On  using  Eqs.(lO)  and  (11)  in  the  wave  function  expression 
(8),  we  can  bring  it  to  the  symmetrized  form 
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In  this  expression  ordering  of  Zj±'^s  is  immaterial.  The  wave  function  (12)  falls  off  exponen¬ 
tially  with  the  distance  \zj^  —  z^_|,  under  the  binding  condition 


m_  >  'm+,  K-  >0;  m_  <  V,  <  0.  (13) 

Clearly,  such  a  wave  function  characterizes  a  bound  state  of  iY  excitons,  as  evident  from 
the  real  exponential  factors  and  an  overall  translation  (the  first  factor),  analogously  to  a 
delocalized  Wannier  excitonic  complex  [6]. 

The  exponential  factors  in  the  bound  state  \^2N  >  can  be  revealed  by  the  dependence 
of  the  intensity-intensity  correlation  function  [3]  on  the  separation  r/  of  two  photon 
counters  detecting  the  field  in  the  structure: 

=  J  <  ^2n\S^{z)S~{z)£'^{z  +  T])S~{z  +  T})\^2N  >  dz,  (14) 

where  the  operator  £  (z)  =  <?^-f(~)  +  is  the  position-dependent  negative-frequency 

field  envelope.  A  laborious  extension  of  the  treatment  in  Ref.  [3]  to  the  two-band  system 
yields  as  a  linear  combination  of  exponential  factors  of  the  form  exp[— (c,-  -t-  c±)|7/|], 

exp[-(c+  +  c_)|7|]  and  exp[-(ci  +  c+/2  -f  c_/2)|7?|],  where  c,  =  2m_m+l^/(m_  - 
c±  =  2|Ka|m±/,^  ,  and  a  =  -P  or  -  according  to  the  Kerr  nonlinearity  sign.  VVe  have 
evaluated  ^^^^(77)  for  iV  =  2  (a  biexciton)  in  a  periodic  structure  consisting  of  Kerr- non  linear 
GaAs  layers  alternating  with  linear  dielectric  layers,  using  the  linear  dispersion  relation  from 
Ref.[9]  and  a  beam  spot  size  p  =  0.5  pra.  The  layer  thicknesses  are  ai(GaAs)  =  0.2338  pm 
and  02  =  0.3044  ^m.  The  refractive  index  of  GaAs  is  rii  =  3.60  and  the  refractive  index 
02  of  linear  dielectric  varies  from  1.00  to  3.574.  GaAs  has  a  nonlinear  susceptibility  yf'  = 
— 2.5482X  10~^°  (cm/V)^  for  light  frequencies  below  the  bandgap  Eg  =  2.1573  x  10^®  s“^  [10]. 
The  potentials  K,  V±  and  effective  masses  m±  have  been  evaluated  for  the  ’’conduction”  and 
’’valence”  bands  bordering  the  second  lowest  band  gap,  for  which  the  binding  condition  in 
Eq.(12)  holds  (as  opposed  to  the  lowest  band  gap  in  this  structure,  for  which  this  condition 
fails).  Figure  1  clearly  manifests  the  non-monotonic  dependence  of  the  multi-exponential  fall- 
off  of  with  Tj  upon  the  refractive  index  ratio  ni/n2.  This  dependence  is  a  unique  feature 
of  the  present  solutions,  whereas  in  quantum  solitons  traveling  through  Kerr-nonlinear  fibers 
[3]  the  exponential  fall-off  of  G^'^^tj)  depends  only  on  the  photon  number  N. 

The  bound-state  energy  eigenvalues  are  found  from  Eqs.(9)  and  (10)  to  have  the  following 
form. 


E2N  =  2Nh{(jjQ  rb  u;2jv), 


(15) 
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where  wq  =  (<^+  +  ‘^-)/2  is  the  center  of  the  band  gap  and  the  upper  (lower)  sign  is  chosen 
according  to  that  of  the  Kerr  nonlinearity.  As  seen  from  Eq.(15),  bound  states  are  associated 
with  discrete  transmission  lines  at  0,2N  =  (Uq  i  cu2W  in  the  band  gap,  whose  progression 
scales  quadratically  with  half  the  photon  number  in  an  incident  pulse  N  =  l2ivl^Ti^2N, 
T  being  the  pulse  duration  and  I2N  the  2iV-photon  component  in  the  power  density.  We 
have  assumed  that  the  incident  light  field  is  in  a  quasi-monochromatic  coherent  state  with 
band-gap  frequency.  Remarkably,  the  discreteness  of  the  transmission  spectrum  will  separate 
out  different  Fock  states  in  the  incident  coherent  state  according  to'  their  relative  weights, 
i.e.,  the  spectrum  will  reflect  the  photon  statistics!  This  discreteness  will  be  revealed  only  if 
|m_  —  is  spectrally  resolvable,  otherwise  the  quasiclassical  limit  of  gap  solitons  will 

prevail.  The  lowest  line  iV  =  1,  which  marks  the  minimal  transmitted  intensity  associated 
with  a  single  exciton,  is  the  farthest  removed  from  the  appropriate  band  edge,  and  the  lines 
are  drawn  nearer  to  the  edge  as  increases,  as  in  the  case  of  ordinary  Wannier  excitons  [6j. 
However,  whereas  Wannier  excitons  have  a  Rydberg  spectrum,  with  an  infinite  number  of 
levels  (which  scale  as  1/fV^),  the  present  bound  spectrum  is  terminated  at  Nmax,  for  which 
Eq.(15)  yields  wo±u;2/v  =  marking  what  can  be  dubbed  the  ’’ionization”  threshold.  Since 
the  wave  vector  K  of  excitons  is  measured  relative  to  2ko,  the  effective  mass  approximation 
requires  K  <  2ko,  so  that  we  can  let  K  =  0  in  Eq.(15). 

The  spectral  features  discussed  above  are  shown  in  Fig.2  for  the  same  structure  as  in 
Fig.l.  Here  too  the  dependence  on  ni/n2  is  strongly  non- monotonic.  A  rather  restricted 
range  of  ni/n2  allows  the  resolution  of  levels  corresponding  to  different  A  on  a  0.1  GHz 
scale.  Using  a  relative  refractive  index  nijn2  —  3.4221  and  a  pulse  duration  r  =  1.0  ps,  we 
have  obtained  the  transmitted  intensities  I2N  =  8-Ff  x  10“'^,  3.07  x  10”^,  7.28  x  10  8.99 

and  27.41  W/cm^,  respectively,  corresponding  to  the  quantum  numbers  A  =  1,  35,  83,  lO'’ 
and  2.66  x  10^  which  is  Nmax-  For  A  >  Nm.ax  the  spectrum  (15)  is  replaced  by  continuous 
bands  of  energy  eigenvalues,  embedded  in  the  allowed  ’’conduction”  and  ’’valence”  bands. 
These  bands  are  analogous  to  the  classical  ” out-gap”  solutions  obtained  above  a  certain 
intensity  threshold  [11]. 

To  sum  up,  we  have  accomplished  the  Bethe-ansatz  solution  for  pairwise  interacting 
quanta  within  the  effective-mass  regime  of  band-gap  propagation  in  ID-periodic  Kerr- 
nonlinear  dielectric  structures.  Our  theory  predicts  a  new  kind  of  collective  excitations  of 
the  electromagnetic  field  dressed  by  such  media,  namely,  optical  multi-exciton  (OME)  com¬ 
plexes  (or  condensates),  which  are  the  quantum  states  associated  with  gap  solitary  waves. 
Their  existence  should  be  manifested  (i)  by  the  discrete  spectrum  of  band-gap  transmission, 
which  reproduces  the  incident  photon  statistics,  i.e.,  spectrally  filters  photon-number  states 
on  a  MHz  scale,  and  (ii)  by  the  multi-exponential  fall-off  of  intensity-intensity  correlations 
on  a  0.1  mm  scale.  Another  distinct  signature  of  OMEs  is  the  strong  non- monotonic  de¬ 
pendence  of  the  above  properties  on  the  periodic  modulation  of  the  refractive  index  in  the 
structure. 

OMEs  should  have  advantageous  stability  properties,  provided  both  linear  and  nonlin¬ 
ear  absorption  are  negligible.  The  stability  properties  of  OMEs  stem  from  their  being  the 
exact  bound  eigenstates  of  the  nonlinear  field  Hamiltonian  which  accounts  for  single-photon 
and  two-photon  processes  in  the  medium.  In  non-absorbing  media,  OMEs  should  be  sta- 
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ble  against  (a)  spontaneous  radiative  recombination,  which  is  the  main  decay  mechanism 
of  electron- hole  excitons  in  semiconductors  [6];  (b)self- focusing  (-defocusing)  by  the  Kerr 
nonlinearity;  (c)pulse  (wavepacket)  dispersion  and  diffraction  in  the  dielectric  structure; 
(d)quantum  fluctuations  in  the  electromagnetic  field.  On  the  other  hand,  in  order  to  ensure 
the  OME  stability  against  thermal  scattering  due  to  thermal  fluctuations  in  the  dielectric 
function,  one  needs  low  temperatures.  If  the  dielectric  structure  is  Raman  inactive,  i.e., 
optical  phonon  scattering  is  absent,  then  thermal  scattering  does  not  affect  the  resolution  of 
OME  lines  on  a  0.1  GHz  scale.  The  ability  to  resolve  OME  lines  corresponding  to  much  lower 
N  would  necessitate  band-structure  design  aimed  at  increasing  the  effective  mass  difference 
|m_  —  m+|  and  the  intraband  nonlinear  potentials  |V±|. 
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FIG.  1.  Dependence  of  the  intensity  correlation  function  G^^\r])  on  the  detector  separation  r), 
for  the  photon  number  2N  =  4. 
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FIG.  2.  Variation  of  the  optical  exciton  frequency  fi2JV  with  the  relative  refractive  index  nxjn^. 
The  insert  displays  the  variation  of  the  band  edges  0;+  and  a;_  with  ni/712  and  cjq  =  (^+  +^-)/2 
is  the  center  of  the  band  gap. 
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We  identify  the  universal  mechanism  responsible  for  superluminal  (faster- 
than-light)  traversal  times  as  well  as  narrowing  of  wavepackets  transmitted 
through  various  non-dissipative  media.  This  mechanism  is  shown  to  be  pre¬ 
dominantly  destructive  interference  between  successive  wavepacket  compo¬ 
nents  traversing  aU  accessible  causally-retarded  paths.  It  strongly  depends  on 
wavepacket  coherence  and  width,  and  can  cause  superluminal  traversal  not 
only  in  evanescent-wave  ’’tunneling”  but  also  in  allowed  propagation. 
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I.  INTRODUCTION 

As  shown  by  a  two-photon  interference  experiment  [1],  a  photon  that  has  tunneled  as  an 
evanescent  wavepacket  through  a  dielectric-mirror  ’’barrier”  appears  to  have  been  delayed 
significantly  less  than  its  ’’twin”  photon  that  has  traversed  the  same  distance  in  vacuum. 
Such  a  delay  has  been  interpreted  as  signifying  ’’superluminal”  (faster-than-light)  barrier- 
traversal  time.  Similar  ’’superluminal”  time  delay  in  tunneling  through  a  dielectric  mirror 
has  now  been  measured  in  a  classical  two- pulse  interference  experiment  [2].  The  latter  ex¬ 
periment  has  also  revealed  a  remarkable  feature,  namely,  that  the  temporal  width  of  the 
transmitted  wavepacket  is  strongly  narrowed  down.  These  intriguing  time-domain  measure¬ 
ments  add  new  insight  to  that  offered  by  earlier  spectral-domain  observations  of  superlumi¬ 
nal  mean-phase  delays,  in  frustrated  total  internal  reflection  [3]  and  waveguide  transmission 
[4,5]  of  electromagnetic  (EM)  evanescent  waves.  Super  luminal  delays  should  also  occur  in 
tunneling  of  massive  particles  through  potential  barriers  [6,7],  which  is  analogous  to  EM 
evanescent- wavepacket  transmission  [8]. 

It  is  always  possible  to  trace  numerically  the  evanescent  wavepacket  evolution  and  com¬ 
pare  its  features  with  different  definitions  of  barrier  traversal  times  [6,7,9]  (see  below).  Nev¬ 
ertheless,  the  mechanism  of  superluminal  time  delays  is  still  obscure  [1]  and  regarded  as  a 
’’poorly  resolved  mystery”  [9].  A  commonly  invoked  notion  is  that  this  mechanism  is  spectral 
reshaping  (filtering)  of  the  transmitted  wavepacket  by  dispersion.  Indeed,  such  reshaping 
explains  pulse  narrowing  and  superluminal  pulse  traversal  in  absorbing  [10]  (or  amplify¬ 
ing  [11])  media,  whose  dispersion  causes  the  faster  spectral  portion  to  be  less  absorbed 
(or  more  amplified)  than  the  slower  one.  Analogous  reshaping  occurs  in  non-relativistic 
electron  wavepackets  which  are  dispersed  in  free  space  before  hitting  the  barrier  [9].  Yet 
why  should  spectral  reshaping  necessarily  yield  superluminal  delays  of  EM  pulses  in  non¬ 
absorbing  structures,  after  propagating  in  (dispersionless)  vacuum?  Is  there  a  common 
mxchanism  for  superluminal  time  delays  and  wavepacket  narrowing,  which  applies  to  both 
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EM  pulses  in  dielectric  structures  and  relativistic  massive  particles  in  potential  barriers 
[4]?  How  is  causality  compatible  with  superluminal  transmission,  particularly  in  the  single¬ 
photon  case  [1]? 

We  purport  to  show  in  this  paper  that  the  above  questions  can  only  be  answered  by 
a  universal  description  of  the  temporal  wavepacket  transmission  as  interference  between 
its  causally-propagating  consecutive  components  [12].  Our  description  reveals,  for  the  first 
time,  the  key  role  of  phase  coherence  in  tunneling,  by  demonstrating  its  dependence  on  the 
coherence  time  (phase  randomization)  of  the  wavepacket.  An  important  corrolary  is  that 
superluminal  time  delays  can  occur  also  in  allowed  propagation,  namely,  propagation  which 
can  only  be  described  by  real  wavevectors  (e.g.,  in  Fabry-Perot  structures)  and  not  only  in 
evanescent- wave  tunneling,  where  complex  wavevectors  can  be  employed  (e.g.,  in  photonic 
band-gap  structures). 


II.  TRANSMISSION  SPECTRUM  AND  MEAN  TRAVERSAL  TIME 


Our  general  framework  assumes  a  classical  EM  pulse  with  field  amplitude  'ipinix.t)  that 
is  normally  incident  from  x  <0  onto  a  dielectric  structure  inO  <  x  <  L,  The  field  amplitude 
'f’trixj  t)  transmitted  through  the  structure  is  measured  at  x  —  L.  It  is  related  to  the  incident 
field  amplitude  at  a;  —  0  [13]  by  convolution  with  the  impulse  response  (X^t) 


roo 

=  /  dr  (j(r)  -  t) 

^  0 


(1) 


where  (7{t)  is  the  Fourier  transform  of  the  spectral  transmission  coefficient  a{ui).  The  defi¬ 
nitions  of  traversal  times  vary  according  to  the  wavepacket  feature  that  is  monitored  [6,7,9], 
We  shall  pick  two  of  them  and  attempt  to  explain  why  they  are  superluminal,  i.e.,  shorter 
than  their  free-space  counterparts  (similar  explanation  of  other  traversal  times  is  deferred 
to  Ref.  14):  (i)  The  mean  traversal  time  is  defined  by 


i-nean  —  ^  {'4’tr'^tr)  (2) 

where  the  overbar  denotes  integration  over  all  times  and  the  pointed  brackets  stand  for  an 
ensemble  average  (required  for  fluctuating  fields).  This  traversal  time  coincides  with  the 
center-of-gravity”  arrival  time  [9]  in  the  case  where  the  propagation  outside  the  barrier  is 
not  dispersive.  We  note  that  tmean  can  be  defined  for  any  shape  of  the  incident  pulse,  (ii) 
The  peak  traversal  time,  tpeafc,  can  be  defined  if  \^ptT\  is  a  smooth  single-peaked  function, 
as  =  0.  Only  in  the  asymptotic  limit  of  a  spectrally  narrow  incident 

pulse,  where  the  stationary  phase  approximation  is  valid,  these  two  traversal  times  coincide 
with  the  so-called  ’phase  time”  [6,7]  tp^ase  —  d(j)Jduj  where  ^(cu)  is  the  phase  of  the  spectral 
transmission  function  5’(cu). 


Ill,  SUPERLUMINAL  TRAVERSAL  TIMES  IN  LAYERED  STRUCTURES 

We  first  show  how  superluminal  effects  occur  in  the  simple  case  of  a  single  dielectric 
layer,  the  region  0  <  x  <  T  is  filled  with  a  (non-dispersive)  medium  of  refractive  index  n2 
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and  embedded  in  an  infinite  medium  of  refractive  index  ni.  The  transmission  coefficient 
o’{uj)  can  be  expanded  in  a  series 


o'(w)  =  (3) 

j=o 

where  each  term  represents  a  causal  path,  corresponding  to  2j  boundary  reflections;  j 
round  trips  through  L  followed  by  transmission  (Fig.  1,  inset).  Here  Tj  =  {2j  +  l)n2Lfc  is 
the  jth  path  traversal  time  and  Cj  =  (1  —  are  determined  by  the  reflection  amplitudes 

A  =  (ni  —  n2)/(ni  +^2)  at  the  boundaries  a;  =  0  and  x  =  L.  Equation  (3)  also  describes  the 
transmission  of  a  Fabry-Perot  etalon  (two  mirrors  separated  by  a  dielectric  medium)  when  A 
is  replaced  by  the  (complex)  reflection  coefficient  of  a  mirror  [14].  By  Fourier-transforming 
the  expansion  (3)  we  obtain  the  impulse  response 

00 

<^(^)  =  (4) 

J=0 

consisting  of  successive  impulses  with  causal  propagation  times  and  decreasing  amplitudes. 

Suppose  that  the  incident  pulse  ■i/j,vi(0,^)  is  of  the  form  ■0m(O,f)  =  5^(0  where 

g{t)  is  a  normalized  Gaussian  of  temporal  width  At  and  the  phase  ({t)  corresponds  to 
an  oscillation  with  frequency  A  and  phase-coherence  time  such  that  = 

gtw(<-<')g-|t-«'|/Tc  tool  will  be  the  following  autocorrelation  function 

r(rj  -  Tfc)  =  Re(V’*n(*5  -  -  n))  = 


=  cos[u;(rj  —  r^)]  exp 


(a  -ny 

8A? 


(5) 


r(rj  —  Tk)  consists  of  the  cosine  of  the  relative  phase  between  the  paths  i  and  j,  weighted 
by  an  exponential  term,  measuring  the  amount  of  overlap  and  phase  correlation  between 
wavepackets  traversing  the  two  paths. 

The  time-integrated  and  ensemble- averaged  transmitted  intensity  ItT  —  {rptr'i’tr)  is  ex¬ 
pressed  in  term  of  this  autocorrelation  function  as 


hr  =11^5  +  11  CjCkT{Tj  -  Tk)  (6) 

If  the  pulse  is  either  very  narrow  (At  C  tq)  or  incoherent  (tc  <C  tq)  then  the  wavepackets 
traveling  along  different  paths  Tj  Tk  are  no  longer  correlated,  due  to  phase  randomization  or 
lack  of  overlap.  The  second  (coherent)  term  in  (6)  is  then  washed  out,  and  the  transmission 
becomes  frequency  independent.  In  the  opposite  limit,  when  the  incident  pulse  is  wide, 
smooth  and  coherent  (transform-limited),  such  that  A(,rc  <C  tq),  strong  interference  takes 
place  between  the  different  overlapping  wavepackets  and  the  second  term  plays  a  major 
role.  When  this  term  is  large  and  negative,  this  means  that  the  interference  is  strong  and 
predominantly  destructive,  due  to  the  negativity  of  the  r(Tj  —  r*,)  terms  with  the  largest 
CjCk  weights.  In  the  single-layer  example,  the  transmission  of  a  wide,  transform-limited  pulse 
becomes  minimal  when  Co{Tj+i  —  tj)  is  an  odd  multiple  of  tt,  yielding  hr  —  [2nxn2l{n\  +nf)]^. 
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The  mean  traversal  time  (Eq.  (2))  through  such  a  structure  is 

:an  ~  ^  CjTj  +  —  ^  CjCk{Tj  +  r^)r(rj  —  Tk)  / 


(7) 


The  first  (diagonal)  term  in  (7),  which  is  predominant  for  an  incoherent  or  temporally 
narrow  incident  pulse,  is  always  larger  than  the  shortest  causal  arrival  time  Tq  =  n2 L /c,  as 
expected  for  Tj  weighted  with  positive  probabilities  {cj/Yljcj).  For  a  wide,  coherent  pulse 
with  carrier  frequency  u)  such  that  the  transmitted  intensity  (Eq.  (6))  becomes  minimal  due 
to  strong,  predominantly  destructive  interference  between  the  paths,  the  second  (coherent) 
term  in  Eq.  (7)  can  become  sufficiently  negative  to  cause  tmean  <  tq.  In  the  present  example 
the  minimal  value  of  tTnean  is  {n-^  T  which  is  less  than  Ljc  if  n2  —  1,  i.e., 

superluminal,  although  we  deal  with  allowed  propagation! 

The  width  of  the  transmitted  pulse  is  given  by  (Af)^  =  -  tmeanY'i^tr'^tr) ■  It  can 

be  shown  that 


=  A?  + 


qj  +  Tk)  r(rj  Tfc)  — 

i,k 


(8) 


In  the  limit  of  total  incoherence,  we  find  that  (Af)^  is  just  the  sum  of  the  squared  widths 
of  the  incident  wavepacket  Aj  and  the  impulse  response  V,- where  c^  = 
1  opposite  limit  of  strong  interference  between  coherent  and  wide  wavepackets 

r(A  ~  Tk)  ~  It  then  follows  from  Eqs.  (5)-(7)  that  Eq.  (8)  becomes 


(At)2  =  A? 
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[ln(/i.(c.))] 


UJ 


(9) 


If  ill  lies  m  a  dip  of  the  spectral  transmission  curve  /jr(cu)  then  d‘^[lnltr]/duj‘^  >  0  and  the 
p^ulse  will  be  narrowed.  The  temporal  narrowing  effect  will  be  most  salient  when  A^  ~ 
aZT  [l^(7tr(w))],  provided  A^  is  large  enough  to  allow  overlap  of  successive  wavepackets.  This 
effect  is  seen  to  be  sensitive  to  coherence:  the  phase  incoherence  which  contributes  only 
to  the  total  spectral  width  of  and  xl>tr  (due  to  the  fluctuating  phase  <f(t)),  exponentially 
diminishing  the  narrowing  in  (8). 

Figure  1  allows  more  insight  into  the  case  of  wide  and  coherent  for  which  ^{r(T,t) 
(thick  curve)  consists  of  overlapping,  destructively  interfering  Gaussians  (thin  curves).  Am¬ 
plitude  suppression  of  the  first  transmitted  wavepacket  |co^,„(t  —  to)|  by  the  next  one  is 
stronger  throughout  its  rear  half  than  throughout  its  forward  half,  since  the  envelopes 
—  Tj)\  with  j  >  1  are  maximal  at  t  >  tq.  Consequently,  the  forward  tail,  which 
has  extended  through  the  structure  already  at  t  <  0,  becomes  the  peak  of  'iptTi  correspond¬ 
ing  to  superluminal  tpgak  and  tmean-  Destruction  of  the  rear  half  of  i/>,„  by  interference  also 
makes  the  transmitted  pulse  narrower,  because  it  consists  mostly  of  the  forward  tail  of 
By  comparison,  an  incoherent  (fluctuating)  input  Gaussian  of  the  same  envelope,  results  in 
a  broad,  intense,  (IV’trP)  with  subluminal  tpeak- 

The  foregoing  results  render  the  essence  of  superluminal  effects  in  the  transmission 
through  any  layered  structure,  since  the  impulse  response  is  then  a  discrete  sum  of  S- 
functions  as  in  (4).  Consider  specifically  the  structure  used  in  Refs.  1  and  2,  which  contains 
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N  periods,  each  consisting  of  two  layers  with  refractive  indices  ni  and  n2,  and  thicknesses 
£i,£2-  We  have  obtair^ed  the  transmission  coefficient  ^^{u})  of  N  such  periods  by  a  recurrence 
relation  from  the  single-period  transmission  and  reflection  coefficients.  Its  Fourier  transform 
can  be  expanded,  as  seen  from  Fig.  2  (inset),  in  the  form 

OO 

=  X]  (10) 

ii  J2-0 


Here  each  coefficient  is  the  sum  of  the  amplitudes  of  all  possible  paths  traversing  the 

Til  and  712  layers  N  +  2ji  and  N  -f  2j2  times,  respectively,  with  causal  delay-times  rji  — 
{N  +  2ji)ni£il2c  {i  =  1,2).  Equations  (6)-(8)  and  the  foregoing  discussion  of  the  single-layer 
(Fig.  1)  fully  apply  to  the  multi-layered  barrier,  on  substituting  Cj  -+  and  Tj  — >  Tji+Tj2. 

We  find  that  (Fig.  2)  an  evanescent  wave  is  a  sum  of  propagating  transmitted  waves  whose 
leading  terms  interfere  destructively  at  band  gap  frequencies.  Correspondingly,  there  is 
constructive  interference  in  the  leading  reflected  waves  (the  Bragg  reflection  condition). 


IV.  SUPERLUMINAL  TRANSMISSION  TIMES  IN  DISPERSIVE  MEDIA 

The  foregoing  results  can  be  extended  to  non-dissipative  dispersive  media  which  are 
characterized  by  a  continuous  impulse  response,  e.g.,  a  dielectric  medium  with  continuously 
varying  refractive  index  or  an  optical  waveguide  [4,9].  The  sums  over  paths  J2j,k  CjCkT{Tj-Tk) 
in  Eqs.  (6)-(8)  should  be  replaced  by  the  integrals  f  dr  /  dr'(T(r)cr(r')r(T  -  r').  Reduced  in¬ 
tensity,  superluminal  time-delay  and  temporal  narrowing  of  pulses  transmitted  through  such 
media  follow,  as  in  layered  structures,  from  destructive  interference  between  components  of 
thtT  at  different  delay  times  r  ^  r',  corresponding  to  predominantly  negative  contributions 
of  (7(r)c7(r')r(T  —  r').  As  an  example,  consider  an  infinite  waveguide,  with  cutoff  frequency 
tUc-  Its  impulse  response  for  transmission  from  x  —  0  to  x  —  L  has  been  found  by  us  to  have 
the  causal  form  [15] 

a{t)  =  6{t  —  Ljc)  —  Uc{Ljc)  [Ji(a;cs)/s]  6{t  —  L/c)  (11) 

where  Ji  is  the  first-order  Bessel  function,  9  is  the  Heaviside  step  function  and  s  = 
If*.  is  coherent  and  temporally  wide  (At,rc  <C  L/c)  we  can  divide  the 
integrals  over  (T(r)cr(r')r(T  —  r')  into  intervals  that  exhibit  strong  cancelations  in  the  ex¬ 
pressions  for  tmearn  hr  and  (Af)^  (the  continuous  limits  of  Eqs.  (6), (7)  and  (8),  respectively) 
provided  that  Co  <  io^.  The  superluminal  delays  observed  in  waveguide  transmission  below 
cutoff  [4]  are  obtainable  from  this  description.  It  is  important  that  massive  relativistic  par¬ 
ticles  in  potential  barriers  obeying  the  Klein- Gordon  equation  fit  the  same  description,  since 
their  energy- momentum  dispersion  is  the  same  as  in  light  in  a  waveguide  [8]! 


V.  CAUSALITY  AND  SIGNAL  TRANSMISSION 

Had  the  peak  of  the  transmitted  wavepacket  carried  any  new  information,  its  arrival 
after  a  superluminal  time-delay  tpeak  would  have  violated  causality.  However,  for  an  analytic 
input  in  Eq.  (1)  the  interfering  forward  tails  of  Tpin{t  —  rj)  which  are  present  at  x  =  T  at 
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t  <  0,  already  contain  all  information  on  the  rest  of  the  pulse  to  follow.  New  information  can 
only  be  transmitted  by  a  non  analytic  (abrupt)  disturbance  of  which  would  be  causally 
delayed.  This  can  be  experimentally  demonstrated  by  switching  off  the  incident  Gaussian 
at  its  peak  (t  =  0)  much  faster  than  (for  ui  in  the  microwave  range,  this  is  achievable 
by  a  sub-picosecond  optical  pulse  that  drastically  changes  the  transmissivity  of  a  dielectric 
medium).  Such  abrupt  switching-off  corresponds  to  ^.•„(a:  =  0,f)  =  e{-t)g{t)  where  6{t)  is 
the  Heaviside  step  function  [16].  The  transmission  of  this  through  any  layered  structure, 
e.g.,  a  single-layer  (Eq.  (4))  yields  by  Eq.  (1) 


Mt)  =  -  t)9{t  -  h  )  (12) 

i 

The  interference  is  unaffected  by  this  switching-off  at  f  <  tq.  Hence,  the  forward  half  of 
looks  the  same  as  for  an  unchopped  Gaussian  and  may  exhibit  a  superluminal  ipea/t-  The 
true  character  of  'thin  is  revealed  only  at  f  =  tq  =  T/c,  when  the  first  transmitted  Gaussian, 
V’m(^  ~  "^o))  vanishes,  causing  'iptr  (thick  curve  in  Fig.  3)  to  drop  below  |cp^,-„(t  —  ri)|.  This 
demonstrates  a  fundamental  point:  The  step-like  decrease  at  successive  tj,  transmits  the 
switching-off  information  in  a  causal  fashion,  whereas  superluminal  features,  such  as  tpeak, 
carry  no  information. 


VI.  PHOTON  DETECTION  AND  SUPERLUMINAL  EFFECTS 

The  discussion  thus  far  has  been  classical,  but  it  can  easily  be  rendered  in  quantal  terms, 
appropriate  for  two-photon  interference  [1]  or  one-photon  detection.  We  must  replace  Tpin  by 
a  field  operator  $,-„^(a:,i)  =  where  are  the  annihilation  operators  of  the 

modes  k  of  the  free  field.  In  the  single-photon  case  the  incident  field  state  is  |1)  =  E«^«a«|0) 
where  g^,  is  a  Gaussian  function  of  a;(«;)  centered  around  Co  and  |0)  is  the  vacuum  state. 
In  this  case  r(r_,'  -  r^)  in  Eq.  (5)  must  be  evaluated  using  (l|^(-)(a;,  t)^(+)(x',  f')|l)  = 
,t')  where  is  the  classical  wavepacket.  Equations  (6)-(8)  are  then  valid,  if 
Itr  is  interpreted  as  the  detection  rate  or  probability  of  a  transmitted  photon.  A  transmitted 
photon  is  likely  to  be  detected  at  ’’superluminal”  times  when  its  Itr  is  peaked  at  ih.Q  forward 
tail  of  the  classical  |z/»,„(T,  t)p,  which  was  already  present  at  the  detector  even  at  t  <  0. 
By  contrast,  a  similar  photon  that  has  propagated  the  same  distance  through  vacuum  is 
characterized  by  the  peak  of  |•!/>^•„(Z,,  i)|^,  which  arrives  latter  (at  t  =  Ljc). 


VIL  CONCLUSIONS 

Our  theory  has  demonstrated,  for  the  first  time,  that  the  universal  mechanism  of  pre¬ 
dominantly  destructive  interference  between  accessible  causal  paths  [12]  is  responsible  for 
transmission  attenuation,  superluminal  delay  times  and  wavepacket  narrowing.  Two  other 
characteristics  of  evanescent  waves  , namely,  exponential  attenuation  and  traversal-length 
independence  of  tmean  or  tps^k  for  opaque  barriers  [2]  can  also  be  explained  in  terms  of  this 
universal  mechanism  [15].  This  theory  overcomes  the  limitations  of  previous  approaches, 
since  it  applies  to  arbitrary  pulse  shapes,  widths  and  coherence  times,  and  explicitly  re¬ 
veals  the  causal  nature  of  their  transmission.  The  understanding  provided  by  this  theory 
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may  open  new  perspectivespn  the  design  of  the  velocity,  intensity  and  shape  of  transmitted 
pulses,  by  manipulating  the  phase  delays  along  the  accessible  paths  in  the  medium. 
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FIG.  1.  A  Gaussian  pulse  transmission  through  an  empty  layer  (112  =  1)  between  two  dielectric 
slabs  with  ni  =  4.  Inset  -  Successively  transmitted  terms  corresponding  to  2j-fold  reflections  in 
the  empty  layer.  Thin  curves  -  successively  transmitted  Gaussians  numbered  by  j  =  0, 1,  2, ...  with 
phase  delays  uTj  —  (2j  +  l)7r.  Dotted  curve:  their  coherent  sum  'iptr  for  narrow  (A^  =  0.03  Ljc). 
Thick  curve:  for  a  wide  i^t  ~  1.5X/c)  exhibits  tmean  ~  '^peak  <  Ljc  (arrow).  Dashed  curve: 

the  square  root  of  the  ensemble- averaged  intensity  when  =  0,5Lfc. 
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FIG.  2.  Transmission  through  5  double-layer  periods  of  length  2i  (ni  =  1.5,  ^2  =  1).  Inset- 
successively  transmitted  terms  due  to  interlayer  reflections.  Impulse  response  spikes  are  shown  at 
t  =  Tji  Tj2^  Solid  curve:  transmitted  Gaussian  envelope  with  carrier  frequency  u)  at  the  center  of 
the  band-gap  and  At  =  1.5T/c  exhibits  t^^ak  ^  0.6i/c. 
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Transmitted  envelope 


FIG.  3.  Same  as  Fig.  1  for  an  incident  Gaussian  pulse  that  is  chopped  at  t  =  0  (Eq.  (12)). 
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